Symmetry, Integrability and Geometry: Methods and Applications 



SIGMA 7 (2011), 056, 41 pages 



Quantum Group U q (sl(2)) Symmetry 

and Explicit Evaluation of the One-Point Functions 

of the Integrable Spin-1 XXZ ChainB 



Tetsuo DEGUCHI and Jun SATO 



Department of Physics, Graduate School of Humanities and Sciences, Ochanomizu University 
2-1-1 Ohtsuka, Bunkyo-ku, Tokyo 112-8610, Japan 



E-mail: deguchi@phys.ocha.ac.jp , jsato @sofia.phys. ocha. ac.jp 



Received October 29, 2010, in final form May 26, 2011; Published online June 10, 2011 



doi:10.3842/SIGMA.2011.056 



Abstract. We show some symmetry relations among the correlation functions of the in- 
tegrable higher-spin XXX and XXZ spin chains, where we explicitly evaluate the multiple 
integrals representing the one-point functions in the spin-1 case. We review the multiple- 
integral representations of correlation functions for the integrable higher-spin XXZ chains 
derived in a region of the massless regime including the anti-ferromagnetic point. Here we 
make use of the gauge transformations between the symmetric and asymmetric i?-matrices, 
which correspond to the principal and homogeneous gradings, respectively, and we send 
the inhomogeneous parameters to the set of complete 2s-strings. We also give a numerical 
support for the analytical expression of the one-point functions in the spin-1 case. 
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1 Introduction 

The correlation functions of the spin- 1/2 XXZ spin chain has attracted much interest during 
the last decades in mathematical physics, and several nontrivial results such as their multiple- 
integral representations have been obtained explicitly [H EJ [3]. The Hamiltonian of the XXZ 
spin chain under the periodic boundary conditions (P.B.C.) is given by 

L 

Kxxz = ^2 ( a f a f+± + a J a J+i + Ao f °f+i) • 



Here <t" (o = X, Y, Z) are the Pauli matrices defined on the jth site and A denotes the anisotropy 
of the exchange coupling. The P.B.C. are given by o^+i = erf for a = X, Y, Z. 

The XXZ Hamiltonian shows the quantum phase transition: the ground state of the XXZ 
spin chain depends on A. For |A| > 1 the low-lying excitation spectrum at the ground state has 
a gap, while for |A| < 1 it has no gap. Here we remark that the quantum phase transition that 
we have discussed is associated with the behavior of the XXZ spin chain in the thermodynamic 
limit: L — > oo. In terms of the q parameter of the quantum group U q (sl2), we express A as 
follows 



*This paper is a contribution to the Proceedings of the International Workshop "Recent Advances in Quantum 
Integrable Systems". The full collection is available at |http:/ /www.emis.de/journals/SIGMA/RA QIS201 0.html| 
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It is often convenient to define parameters rj and ( by q = expr/ with r\ = i£. Here we have 
A = cosh?] = cos£. In the massive regime A > 1, we set rj > 0. In the massless regime 
— 1 < A < 1, we set r? = i£ where £ satisfies < £ < tt. Here, the XXX limit is given by 
i] — > +0 or ( — > +0. Here we remark that the XXZ Hamiltonian can be derived from the 
i?-matrix of the affine quantum group with q parameter, U q {sl2)'- we derive the .R-matrix by 
solving the intertwining relations, construct the XXZ transfer matrix from the product of the R 
matrices, and then we derive the XXZ Hamiltonian by taking the logarithmic derivative of the 
XXZ transfer matrix. Thus, the q parameter of the affine quantum group is related to the 
ground state of the XXZ spin chain through A. 

The multiple-integral representations of the XXZ correlation functions were derived for the 
first time by making use of the q- vertex operators through the affine quantum-group symmetry 
in the massive regime for the infinite lattice at zero temperature [H[2]. In the massless regime 
they were derived by solving the q-KZ equations [51 [6]. Making use of the algebraic Bethe- 
ansatz techniques [TJ [lj [HJ [10], the multiple- integral representations were derived for the 
spin-1/2 XXZ correlation functions under a non-zero magnetic field [TT]. Here, they are derived 
through the thermodynamic limit after calculating the scalar product for a finite chain. The 
multiple- integral representations were extended into those at finite temperatures [12], and even 
for a large finite chain [13] . Interestingly, they are factorized in terms of single integrals [14j . We 
should remark that the multiple-integral representations of the dynamical correlation functions 
were also obtained under finite-temperatures |15J . Furthermore, the asymptotic expansion of 
a correlation function of the XXZ model has been systematically discussed [16]. Thus, the 
exact study of the XXZ correlation functions should play an important role not only in the 
mathematical physics of integrable models but also in many areas of theoretical physics. 

Recently, the form factors of the integrable higher-spin XXX spin chains and the multiple- 
integral representations of correlation functions for the integrable higher-spin XXX and XXZ 
chains have been derived by the algebraic Bethe-ansatz method [171 [T8l \19\ l20l [2T] (see also [22] ) . 
The spin-1 XXZ Hamiltonian under the P.B.C. is given by the following [23] : 

Ns f - 1 
^spin-i xxz = J/ASj • Sj+i ~ (Sj ■ S j+1 ) 2 - -(q-q l ) 2 [SjS z j+l - (5|5| +1 ) 2 + 2(S z j ) 2 ) 

i=i 1 

— (q + q 1 — 2)[(SjSj+i + Sj +1 )S z S z +1 + + ^j^+i)!- 

Furthermore, the multiple-integral representations have been obtained for the correlation func- 
tions at finite temperature of the integrable spin-1 XXX chain [21]. The solvable higher-spin 
generalizations of the XXX and XXZ spin chains have been derived by the fusion method in 
several references [251 E3 EH [281 121 EH ED E2] . In the region: < ( < vr/2s, the spin-s ground- 
state should be given by a set of string solutions [33|. 134] . Furthermore, the critical behavior 
should be given by the SU(2) WZWN model of level k = 2s with central charge c = 3s/(s + 1) 
[HgElESEgEgililHn For the integrable higher-spin XXZ spin 

chain correlation functions have been discussed in the massive regime by the method of (/-vertex 
operators [Ml |49l EUJ EHJ EU |52] . 

The purpose of this paper is to show some symmetry relations among the correlation functions 
of the integrable spin-s XXZ spin chain by explicitly calculating the multiple-integral representa- 
tions for the spin-1 one-point functions. Associated with the quantum group U q (sl{2)) symmetry, 
there are several relations among the expectation values of products of the matrix elements of 
the monodromy matrices. For the spin-1 case, we confirm some of them by evaluating the 
multiple integrals analytically and explicitly. Here we should remark that the derivation of the 
multiple-integral representations for the spin-s XXZ correlation functions given in the previous 
papers [TH [201 EI] was not completely correct: the application of the formulas of the quantum 
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inverse-scattering problem was not valid |53|, [5l] . We thus review the revised derivation [53\ [53] 
in the paper. The spin-s correlation function of an arbitrary entry is now expressed in terms of 
a sum of multiple integrals, not as a single multiple integral. Furthermore, we show numerical 
results which confirm the analytical expressions of the spin-1 one-point functions. 

Let us express by (E 00 ), (E n ) and (E 22 ), the expectation values of Sf = 1, Sf = and 
Sf = — 1, respectively, where Sf denotes the Z-component of the spin operator defined on the 
first site. Then, we have the following: 



(E 22 ) = {E 



00 \ 



C — sin ( cos £ 
2Csin 2 C ' 



(E 



cos £(sin C — C cos C) 
Csin 2 C ' 



We shall show the derivation of (E 00 ), (E 11 ) and (E 22 ), in detail. Here we remark that the 
expressions of (E 22 ), the emptiness formation probability, and have been reported in [20] 

without an explicit derivation. In fact, although the derivation was not completely correct, the 
expressions of the spin-1 one-point functions are correct [531 154] , Here, the quantum group 
symmetry as well as the spin inversion symmetry play an important role, as we shall show 
explicitly in the present paper. 

It is nontrivial to evaluate the multiple integral representations of the XXX and XXZ models 
analytically or even numerically. Let us now return to the spin- 1/2 case. Boos and Korepin 
have analytically evaluated the emptiness formation probability P{n) of the XXX spin chain for 
up to n = 4 successive lattice sites [55]. Performing explicit evaluation of the multiple integrals, 
they successfully reproduced Takahashi's result that was obtained through the one-dimensional 
Hubbard model [56]. The method was applied to all the density matrix elements for up to 
n = 4 successive lattice sites in the XXX chain |57j and also in the XXZ chain [581 [59"1 160] . 
Furthermore, the algebraic method to obtain the correlation functions of the XXX chain based 
on the qKZ equation has been developed [61] and the two-point functions up to n = 8 have been 
obtained so far [62], [631 E31 [65] • At the special anisotropy A = 1/2, some further results have 
been shown for the correlation functions through explicit evaluation [661 WJ\ 1681 [69] . 

The paper consists of the following. In Section [2] we review the Hermitian elementary matri- 
ces [20], and give the basis vectors and their conjugate vectors in the spin-1 case as an illustrative 
example. We also show a formula for expressing higher-spin local operators in terms of spin-1/2 
local operators in the spin-1 case, which plays a central role in the revised method [53[ 154] . In 
Section [3] we summarize the notation of the fusion transfer matrices and the quantum inverse 
scattering problem for the spin-s operators. For an illustration, in Section [4[ we show some 
relations among the expectation values of the Hermitian elementary matrices in the spin-1 XXX 
case and then in the spin-1 XXZ case. In particular, we show the spin inversion symmetry. We 
also show the transformation which maps the basis vectors of the spin-1 representation 
constructed in the tensor product of the spin-1/2 representations <S> V^ 1 ' to the basis of the 
three-dimensional vector space C 3 . The former basis is related to the fusion method, while the 
spin-1 XXZ Hamiltonian ([l.ip is formulated in terms of the latter basis. In Section [5j we review 
the revised multiple-integral representations of correlation functions for the integrable spin-s 
XXZ spin chain |53[ [54] . Here we remark that necessary corrections to the previous papers [19] 
and [20] are listed in references [20] and [21] of the paper [54], respectively. In Section El we ex- 
plicitly calculate the multiple integrals of the one-point functions for the spin-1 XXZ spin chain 
for a region in the massless regime. We show some details of the calculation such as shifting the 
integral paths. In Section [7] we show that the numerical estimates of the spin-1 one-point func- 
tions obtained through exact diagonalization of the spin-1 XXZ Hamiltonian (jl.ip are consistent 
with the analytical expressions of the spin-1 one-point functions. Thus, we shall conclude that 
the analytical result of the spin-1 one-point functions should be valid. 
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2 The quantum group invariance 

We construct the basis vectors of the finite-dimensional spin-^/2 representation of the quantum 
group U q (sl2) in the tensor product space of the spin-1/2 representations, and introduce their 
conjugate vectors. In terms of the basis and conjugate basis vectors we formulate the spin-£/2 
elementary matrices which have only one nonzero element 1 with respect to entries of the basis 
and conjugate basis vectors. We then illustrate an important formula for reducing the spin-^/2 
elementary matrices into a sum of products of the spin-1/2 elementary operators. 



2.1 Quantum group U q (sl 2 ) 

Let us introduce the quantum group Uq^slz) in order to formulate not only the /2-matrix of 
the integrable spin-s XXZ spin chain algebraically but also the higher-spin elementary matrices, 
by which we introduce correlation functions. Here we remark that the correlation functions of 
the spin-s XXZ spin chains are given by the expectation values of products of the higher-spin 
elementary matrices at zero temperature. 

The quantum algebra U q {sl2) is an associative algebra over C generated by X^, with 
the following relations [7011711172"] : 

KK~ l = K~ X K = 1, KX ± K~ l = q ±2 X ± , \X + ,X~] = K ~ K , 

q - q- 1 

The algebra U q (sl2) is also a Hopf algebra over C with comultiplication 

A(x + ) = x + ®i + k ® x + , A(r) = r®r 1 + i(x)r 1 a(k) = k®k, 

and antipode: S(K) = K~\ S{X + ) = -K~ 1 X + , S(X~) = -X~K, and coproduct: e(X ± ) = 
and e(K) = 1. 



2.2 Basis vectors of spin-£/2 representation of t/ 9 (s/ 2 ) 

We introduce the g-integer for an integer n by [n] q = (q n — q~ n )/(q — q -1 )- We define the 
(/-factorial [ra]J for integers n by 

[n] q \ = [n] q [n-l] q ---[l] q . 

For integers m and n satisfying m > n > we define the (/-binomial coefficients as follows 



m 



[m - n] q \[n] q \' 



Let us denote by \a) for a = 0, 1, the basis vectors of the spin-1/2 representation Here 
we remark that and 1 correspond to f and \,, respectively. In the Ith. tensor product space 
(yi 1 )^®? we construct the basis vectors of the (I + l)-dimensional irreducible representation of 
U q (sl2), | \£, n) for n = 0, 1, . . . , £, as follows. We define the highest weight vector | \£, 0) by 

\\£,0) = |0)i ® |0) 2 ®---® |0> € . 

Here \a)j for a = 0, 1, denote the basis vectors of the spin-1/2 representation defined on the 
jth position in the tensor product (T/W)®^. We define \\£,n) for n > 1 and evaluate them as 
follows \M 

n+«2H Hn—ni+n(n—l)/2 



\\£,n) = (A^Hx-)) n \\£,0)j^ = £ «V"^|0>9 

L ig ' i<h<-<i n <e 
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Here denotes the Pauli spin operator a acting on the jth component of the tensor product 
(yW)®^: we have a~ = (g> a - (g "W e define the conjugate vectors explicitly by 

the following: 

(An|| = 



? ra(*-n) ^ __ a + q i 1 +-+i n -ne+n{n-l)/2_ 

i i<h<-<i n <e 

It is easy to show the normalization conditions [19]: (£,n\\ \\£,n) = 1. Let us define F(£,n) by 



g-n(/-n)_ 

We have (\\£,n)) t \\£,n) = F(£,n), and hence {£,n\\ = (\\£, ri))* /F(£, n). Here the superscript t 
denotes the matrix transposition. 

In the massive regime where q = exp 77 with real rj, conjugate vectors {£, n\ | are also Hermitian 
conjugate to vectors \\£, n). 

2.3 Affine quantum group Uq^sl^) 

In order to define the i2-matrix in terms of algebraic relations we now introduce the affine 
quantum group U q (sl2)- It is an infinite-dimensional algebra generalizing the quantum group 
U q {sl 2 ). _ ' 

The algebra Uq^sfa) is an associative algebra over C generated by Xf ,Kf for i = 0,1 with 
the following defining relations: 

K i K~ 1 = Kr 1 K i = l, K.XfKT^ = q ±2 X±, K.XfKr 1 = q ^xf i^j, 

(x±fxf - [^{xffxfxf + [^xfxfixff - Xf(Xff = 0, i + j. 

The algebra U q (sl2) is also a Hopf algebra over C with comultiplication: 

A(X+) = X? ® 1 + Ki ® Xf , A(Xr ) = X~ ® ifr 1 + 1 g, ir , A(i^) = <g iQ, 

and antipode: S(Ki) = Kr 1 , SpQ) = -K^X+, S(Xr) = -XrK,, and counit: e{Xf) = 
and e{Ki) = 1 for i = 0, 1. 

The quantum group U q (sl2) gives a Hopf subalgebra of U q {sl2) generated by X^, Ki with 
either i = or i = 1. Thus, the affine quantum group generalizes the quantum group U q {sl2)- 

2.4 Evaluation representations with principal and homogeneous gradings 

We shall introduce two types of representations of U q {sl2)'- evaluation representations associa- 
ted with principal grading and that with homogeneous grading. The former is related to the 
symmetric i?-matrix which leads to the most concise expression of the integrable quantum spin 
Hamiltonian, while the latter is related to the asymmetric i?-matrix R + {u) which we shall define 
in Section T3.2I and suitable for an explicit construction of representations of the quantum group. 
Here and hereafter we denote by X^ and K the generators of U q (sl2)- 

Let us now introduce a representation of U q {sl2) associated with homogeneous grading [2]. 
With a nonzero complex number A we define a homomorphism of algebras <p^'- U q {sl2) — > 
U q (sl2), as follows. 

<P?V$) = e ±A * T , <^W) = e ±X X±, $\K ) = K~\ <pf{K x ) = K (2.1) 
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Thus, from a given finite-dimensional representation (jr^\V^) of the quantum group U q {sl2), 
we derive a representation of the affine quantum group U q {sl 2 ) by n^(ip^\a)) for a G U q (sl 2 ), 

where tpf '(•) is given by (|2.ip . We call it an evaluation representation of the affine quantum 
group; more specifically, the spin-^/2 evaluation representation with evaluation parameter A 
associated with principal grading. We denote it by (tt^ , V^(X)) or V^ tp \X). 

Similarly in the case of principal grading, we now introduce a representation associated with 
homogeneous grading [2]. With a nonzero complex number A we define a homomorphism of 

algebras <fx'- U q {sl 2 ) — > U q {sl2) by the following: 

riV±) = e ±2X XT, ^ +) (Af ) = X± 4 +) (K ) = K~\ ^(K,) = if. (2.2) 

From a given finite-dimensional representation (ir^ ,V^) of the quantum group [7<j(sZ 2 ) we 
derive a representation of the affine quantum group U q (sh) by TT^\ip^\a)) for a G U q (sl 2 ), 
where ip^~\-) is given by (|2.2p . We call it the spin-^/2 evaluation representation with evaluation 
parameter A associated with homogeneous grading. We denote it by (vr^ +) , V^(A)) or V^ + \X). 

2.5 Defining relations of the i?-matrix 

Let us now define the i?-matrix for any given pair of finite-dimensional representations of the 
affine quantum group U q (sl 2 )- Let (tti,Vi) and (tt 2 ,V 2 ) be finite-dimensional representations 
of U q (sl 2 ). We define the i?-matrix R\ 2 for the tensor product Vi<8> V2 by the following relations: 

7Ti ® 7T 2 (r o A(a)) #12 = i?i 2 vri ® vr 2 (A(o)) , a£U q (sT 2 ). (2.3) 

Here r denotes the permutation operator: r(a ® 6) = b® a for a, 6 G U q {sl 2 ). 

For an illustration, let us write down relations (|2 . 3[) of the i?-matrices associated with eva- 
luation representations. We call them intertwining relations. Associated with principal grading 
we have for o = X^ , Xf 1 and ifi, respectively, the following relations: 

R^(Xi - X 2 )(e Xl X- ® 1 + e x% K~ x ® X") = (e A2 l ® X" + e Al X" ® if-^i^Ai - A 2 ), 
4?( A i " A 2 )(e" Al X+ ® if + e^l ® X+) = (e~ A2 if ® X + + e~ Xl X + ® (Ai - A 2 ), 

#i?(Ai - A 2 )(e Al X+ ® 1 + e Aa if ® X+) = (e A2 l ® X+ + e Al X+ ® if)i?g(Ai - A 2 ), 
i?S§(Ai - A 2 )(e~ Al X- » if" 1 + e~ A2 l ® X") 

= (e^if" 1 ® X~ + e~ Xl X~ l)R { i 2 ] (Xi - A 2 ), 
R ( ^(X 1 - X 2 )K ®K = K® KB${\ X - A 2 ). (2.4) 

Associated with homogeneous grading we have 

R{V(^i - A 2 )(e 2Al A~ ® 1 + e^if" 1 ® X ) 

= (e 2Aa l® X- + e 2Xl X~ ® K'^R^iXi - A 2 ), 
R$ (Ai - A 2 ) (e- 2Xl X+ ® if + e~ 2X2 1 ® X+) 

= (e~ 2X2 K®X + + e" 2Al X+ ® l)^^! " A 2), 
4t } (Ai -X 2 ){X + ® 1 + if ® A+) = (1® A + + A+®if)4 2 ) (Ai- A 2 ), 

^(Ai-AsKr^r^i^r) = (if- 1 ® x~ + x~ \)R^ 2 \x x - a 2 ), 

i?^ } (Ai - A 2 )if ® if = if ® if ^(Ai - A 2 ). (2.5) 
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Here Ai and A2 correspond to the "string centers" of the sets of the evaluation parameters 
associated with the evaluation representations n\ and tt2- We have Ai = £1 — {£ — l)r//2, if 7Ti is 
given by the spin-£/2 evaluation representation derived from the tensor product (T/W)®^ with 
complete ^-string Wj for j = 1,2, . . . ,£. Here we shall define complete strings in Section 13.61 

We can show that the solution of intertwining relations (|2.3p is unique. We may therefore 
define the i?-matrix in terms of relations (12.3f) . 

We remark that relations (|2.4p for the evaluation representation associated with principal 
grading are mapped into those of (|2,5p associated with homogeneous grading through a similarity 
transformation, which we call the gauge transformation. We shall formulate it in Section 13.41 



2.6 Conjugate vectors and Hermitian elementary matrices 

In order to construct Hermitian elementary matrices in the massless regime where q is complex 
and j q\ = 1, we now introduce another set of dual basis vectors [20J. For a given nonzero integer £ 

we define {£, n\ | for n = 0, 1, . . . , n, by 
-1 



n 



E (<>K~< 



+ . . . n + n -(h+-+in)+n£-n(n-l)/2 



l<i 1 <---<i n <£ 



They are conjugate to \\£, n): (£,m\\\\£,n) = 5 m ^ n . Here we have denoted the binomial coeffi- 
cients for integers £ and n with < n < £ as follows 

A = £\ 
nj (l — n)\nV 

We now introduce vectors \\£, n) which are Hermitian conjugate to {£, n\\ when \q\ = 1 for 
positive integers £ with n = 0,1,... ,£. Setting the norm of \\£,n) such that (£, n\\\\£, n) = 1, 
vectors \\£,n) are given by 

||4 n) = Yl a h ' " arJO)q~^+-+ i «)+ ni - n ( n - i y 2 



l<ii<-<i n <£ 

We have the following normalization conditions: 

-2 



n 



-n(e-n) 



11 



(£,n\\\\£,n) 



'£ 




n 









for n = 0, 1, ... , £. 



In the massless regime where q is complex with \q\ = 1, we define elementary matrices 
£m,n(£+) by 

E m ^ e +) = \\£^m)(£,n\\ for m,n = 0,1, . . . ,£. 

In the massless regime matrix \\i, n)(£,n\\ is Hermitian: (\\£,n}(£,n\\^ = \\£,n)(£,n\\. How- 
ever, in order to define projection operators P such that PP = P, we have formulated vectors 

Associated with principal grading we define the spin-£/2 symmetric elementary matrices 
by pi EI] 



&<Md = \\l,i)(l t 3\ 



' F(£,j) 
F(£,i) 



for i,j = 0, 1, 
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2.7 Projection operators 

We define the projection operator acting on from the 1st to the ith. tensor-product spaces by 

I 



P < S..4 = Y,\\^ n )^ n \Y (2-6) 

?1=0 



~(i) 

We introduce another projection operator P 12 ,,^ as follows 



P[l. t = Y,\\^ n )^ n \V (2-7) 

n=0 

The projector P^...* i s idempotent: (P^...^) 2 = Pi2---r I n the massless regime where q is complex 
with \q\ = 1, it is Hermitian: (P^) 1 " = P^-r ^om (g^D and (1277]) . we show the following 
properties: 

pW pW _ pW fo OA 

r 12--^ r 12--^ — - r 12---£' 

pit) pM _ pM (2 n\ 

2.8 Spin-s elementary matrices in terms of the spin-1/2 elementary matrices 

Let us denote by e a,b such 2-by-2 matrices that have only one nonzero matrix element 1 at the 
entry (a, b) for a, b = 0, 1. We call them the spin-1/2 elementary matrices. We denote by e°j ,b 
the elementary matrices e a ' b acting on the jth component of the tensor product (lA 1 ))®^. 

Let us introduce variables e' a and which take only two values and 1 for a,/3 = 1, 2, . . . , i. 
We define diagonal two-by-two matrices $j by &j = diag(l, exp(wj)) acting on vj 1 ^ for j = 
0,1, ... ,L. Here u>j (1 < j < L) are called the inhomogeneous parameters of the spin-1/2 XXZ 
spin chain, and we set wo = Ao (see also Section [3/3]) . We define the gauge transformation 

by a similarity transformation with respect to the matrix xoi—L = ^o^i • • • 3>l- Here, we put 

(£) 

inhomogeneous parameters Wj with the complete ^-strings such as w^k_i^ + j = w^ k __^ + - = 
£fc ~~ (j ~~ I) 7 ? f° r j = 1) 2, • • • , I and k = 1, 2, . . . , N s . Then, we can show the following relation. 

Proposition 1 ([53, 54J). The spin-£/2 symmetric elementary matrices associated with principal 
grading are decomposed into a sum of products of the spin-1/2 elementary matrices as follows 

1/2 

e -(H)(6-(M),/2)^ £ Xl2 ..,e^ • • • ef E %l. t (2.10) 

Here the sum is taken over all sets of eps such that the number of integers /3 satisfying 8/3 = 1 
for 1 < (3 < £ is equal to j. We take a set of e' a s such that the number of integers a satisfying 
e' Q = 1 for 1 < a < i is equal to i. The expression (I2.10P is independent of the order of e' a s with 
respect to a. 

The formula f|2. 10|) plays a central role in the revised derivation of the spin-£/2 form factors 
and the spin-£/2 XXZ correlation functions |53[ [54]. We shall derive (I2.10P in Appendix lAl 
We recall that the derivation of the multiple-integral representations of the integrable spin-s 
XXZ spin chain given in the previous papers [191 I2Q1 [21] was not completely correct [531 El] ■ 
In fact, the transfer matrix becomes non-regular at A = and hence the straightforward 

application of the QISP formula was not valid. 



jff,j{l.p) — ( 


~i 




'£ 






i 




j_ 
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2.9 Example: spin-1 case 

We shall show reduction formula (|2.10p for the spin-1 case. 
The spin-1 basis vectors ||2,ra) (n = 0, 1,2) are given by pjJ 



(2,0|| = (+ + 



[2], 



l|2,0) = |++), 

||2,l) = | + -)+g- 1 |-+), (2,1| 
l|2,2) = |-->, (2,2|| = (- 

Here | H — ) denotes |0)i ® 1 1)2 , briefly. The conjugate vectors ||2, n) (n = 0, 1,2) are given by 



|2,0) = |++), 



(2,0|| = (+ + 



[2], 



\2,l) = (\+-)+q\-+))^, 

M = !-->. 



(2,1| 

(2) 



Let us derive the projection operator P 12 . Explicitly we have 



|2,1)(2,1| 



(I + ")(+- I + ")(- + I +3I-+X+-I + I-+X- + I) 



1 / 0.0 1,1 . 
2^1 e 2 + 



-1 0,1 1,0 1 1,U U,l . 1,1 u,u\ 

e{ e 2 ' -H^' e 2 ' + e 2 ' j. 



(2.11) 



Here we remark that in the massless regime where q is complex with \q\ = 1, operator ||2, 1)(2, 
is Hermitian while ||2, 1)(2, 1|| is not. As a four-by-four matrix we express P 12 by 

/ 1 \ 

1/2 g-72 

1 g/2 1/2 



5(2) 
12 



|2,0)(2,0|| + ||2,1)(2,1|| + ||2,2)(2,2|| 



\ 







(2.12) 



1 / 



[1,2 



Here the symbol [1, 2] at the bottom of the 4x4 matrix of (12.12P denotes that the matrix acts on 
the tensor product space (g) . We note that operator \-\ — )( — h | corresponds to e^e^' 
in (|2.1ip . which gives the entry of (1,2) in the four-by-four matrix of (j2. 12|) ; i.e., the element in 
the 2nd row and 3rd column. 

For an illustration, let us show reduction formula (|2.10p for the spin-1 case. With e' x = and 
e' 2 = 1, reduction formula (I2.10p for i = j = 1 reads 

EVW) = pX) {2 , 1|| = pV X n(e°A 1 / + e^e^Xll (2-13) 
First, it is straightforward to show 



0,0 1,1 -1 
Xi2e a e 2 X12 



0,0 1,1 



0,1 1,0 -1 
XviH e 2 X12 



-1 0,1 1,0 



-l c-2 ' A12 c l c-2 A12 — V c l c 2 

Then, in terms of the four-by-four matrix notation we have 

/ \ 

1 q- 1 




0,0 1,1 



+ 



1 0,1 1,0 
e l e 2 



(2.14) 



\ / 



[1,2] 



Here e^e^' 1 corresponds to the element in the 2nd row and 2nd column of the 4x4 matrix (|2.14p . 
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Multiplying f)2 . 12|) by (|2.14p and making use of (|2,lip . we have the following relation: 

P«(^ 1 + g- 1 ef 1 4'°) = M(2 > l||. 
We have thus confirmed reduction formula ()2.10j) for £ = 2 and i\ = ji = 1, as shown in (|2.13p . 

3 Fusion transfer matrices and higher-spin expectation values 

We construct the monodromy matrices of the integrable higher-spin XXZ spin chains through the 
fusion method. We then evaluate the form factor of a given product of the higher-spin operators 
by reducing them into a sum of products of the spin-1/2 operators and calculate their scalar 
products of the spin-1/2 operators through Slavnov's formula. When we reduce the higher-spin 
operators, we make use of the fusion construction where all the elements are constructed from 
a sum of products of the spin-1/2 operators multiplied by the projection operators. 

3.1 Tensor product notation 

(2s) 

Let s be an integer or a half-integer. We shall mainly consider the tensor product V± ® • • • ® 
vjq^ of (2s + l)-dimensional vector spaces vj with L = 2sN s . Here Vj^ have spectral 
parameters Xj for j = 1, 2, . . . , N s . We denote by E a,b a unit matrix that has only one nonzero 
element equal to 1 at entry (a, b) where a, b = 0, 1, . . . , 2s. For a given set of matrix elements 
Afr'g for a, b = 0, 1, ... , 2s and a, = 0, 1, . . . , 2s, we define operators Aj & for 1 < j < k < N s 
by' 

a,b=l a, {) 

® Ef ® lf°\ ® • • • ® 1^1 ® Ef ® /gfj ® • • • ® 4 2s ). (3.1) 
In the tensor product space, (v( 2s ^)® Ns , we define _g™> n ( 2su ') f or i = 1 ; 2, . . . , iV s and u; = +,jp 

by 

gm,n(2sw) _ ^ j(2s)^ ^ gm.,n(2sw) ^ ^j(2a)\®{Ns-i) 

The elementary matrices E n,n ^ 2sw ^ for n = 0,1,..., 2s and to = +,p, are Hermitian in the 
massless regime. 

3.2 Asymmetric and symmetric i?-matrices 

Let us introduce the i?-matrix of the XXZ spin chain [TJ El EJ [TT]. Let Vi and V2 be two- 
dimensional vector spaces. We define the .R-matrix Rf 2 acting on V\ ® V2 by 

R+(\ 1 -\ 2 ) = R + (v)? d er®<% 

a,b,c,d=0,l 

where u = \\ — A2, b(u) = sinhu/ sinh(ii + rj) and ^{u) = exp(±u) sinh^/ sinh(u + rj). 

We remark that the R + {\\ — A2) is compatible with the homogeneous grading of U q (sl2). In 
fact, it is straightforward to see that the asymmetric i?-matrix satisfies the intertwining relations 



/ 1 

b{u) 

c+{u) 

\ 





b{u) 
1 / 





c~(u) 



(3.2) 



[1,2] 
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associated with homogeneous grading (|2.5p for the tensor product of the spin- 1/2 representations 

Of U q (sh), V^^V}^. 

We denote by R&^u) or simply by R(u) the symmetric /2-matrix where c^{u) of (|3.2p are 
replaced by c(u) = sinhr// sinh(it + rj) |19j . The symmetric i?-matrix is compatible with evalua- 
tion representation associated with principal grading for the affine quantum group Uq^slz) [19] . 
Hereafter we express R + and Ry?' by R( lw > with w = + and p, respectively. 

3.3 Monodromy matrix of type (1, 

We now consider the (L + l)th tensor product of the spin-1/2 representations, which consists of 
the tensor product of auxiliary space Vq and the Lth tensor product of quantum spaces vj 

for j = 1, 2, . . . , L, i.e. F (1) ® (v}^ ® ■ ■ ■ <g) V^ 1} ) . We call it the tensor product of type (1, 1® L ) 
and denote it by the following symbol: 



(1,1®*) = (1,1,1,...,!). 

Applying definition (|3.ip for matrix elements R(u)^ of a given i?-matrix such as R( lw ^ with 
u> = + and p, we define i?-matrices Rjk(\j, A&) = Rjk(Xj — A&) for integers j and A; with 
< j < k < L. For integers j, A; and £ with 0<j<k<£<L, the i?-matrices satisfy the 
Yang-Baxter equations 

Rjki^j — Afc)i?^(Aj - \e)Rke{\k - ^e) = Rkii^k — ^e)Rje(^j — ^e)Rjk{^j — A^). 
We define the monodromy matrix of type (1, l® i ) associated with homogeneous grading by 

t o,i2-"l( a o; wi,w 2 ,. ■ -,w L ) = Rq L (Xo - w L ) ■ • -i?o" 2 (Ao - w 2 )Rq X {\q - w x ). 

Here we have set Xj = Wj for j = 1,2, . . . , L, where Wj are arbitrary parameters. We call 
them inhomogeneous parameters. We have expressed the symbol of type (1, 1® L ) as (1, 1) in 
superscript. The symbol (1,1+) denotes that it is consistent with homogeneous grading. We 
express operator-valued matrix elements of the monodromy matrix as follows 



C^ L {\-{ Wj } L ) D^ L (X;{ Wj } L ) 



Here {w^l denotes the set of L parameters, w\, w 2 , ■ ■ ■ , wl. We also denote the matrix elements 
of the monodromy matrix by [^qi^^A; {wj}i)] a ,fe for a, b = 0, 1. 

3.4 Gauge transformations 

We derive the monodromy matrix consistent with principal grading, ^o^.^A; {wj}z), from 



that of homogeneous grading via similarity transformation Xoi- L as follows [19] 

= ( X v*..LA^. L {K{w j } L )x£.. L e-^xi2-LB^l L (\;{w 3 } L )x x l. L \ 
\ e x o Xl2 ... L cg p) L (X;{w j } L )xu 1 .. L X12-lD^ l {X; {wj} L )xn-L J ' 

Here we recall that Xoi—l = ^o^i " " " an( i $j are given by diagonal two-by-two matrices 
<3?j = diag(l, exp(wj)) acting on vj for j = 0, 1, . . . , L, and we set wq = Ao- In [T9] operator 
7l (1+) (A) has been written as ^4 + (A). 



12 



T. Deguchi and J. Sato 



We now introduce the gauge transformation for the spin-s representation [5l]. We define 
diagonal matrix <&( 2s \w) on the basis vectors ||2s, n) as follows: 

&( 2s \w)\\2s,n) = exp(nw)\\2s,n) for n = 0, 1, . . . , 2s. 

We denote by Q^ s \w) the matrix & 2s \w) defined on the jth. component of the tensor product 
vj 2 ^ ® • • • (g) vjj . We define Xi%..n ac thig on the quantum space V} 2s) <g> • • • <8> Vf? a s) by 

(2s) ^(2s)/ A \ ^(2s)/ A \ 

We express A;, as = ^ — (2s — 1)77/2 for 6 = 1, 2, ... , iV s . Here denote the inhomogeneous 
parameters of the spin-s XXZ spin chains, which will be given in equation (|3.4p of Section [3^6l We 
note that A;, corresponds to the string center of the 2s-string, — (/3 — 1)77 with /3 = 1, 2, . . . , 2s, 
for each 6 satisfying 1 < b < N s . 

3.5 Projection operators through fusion 

Let V\ and V 2 be the (2s + l)-dimensional vector spaces. We define permutation operator IIx 2 
by 

Ul : 2Vl (8> V2 = V 2 ® Vl, Vl e Vl, D2 G V2. 

In the case of spin-1/2 representations, we define operator Rf 2 (X\ — A2) by 
R+ 2 (X 1 -X 2 ) = U lt2 Rf 2 (X 1 -X 2 ). 

We now introduce projection operators P 12 .« f° r £ > 2. We define P 12 by P 12 = Rf 2 (il)- 



For £ > 2 we define projection operators inductively with respect to £ as follows [TIJ 132] 

The projection operator P 12 ...£ gi ves a (/-analogue of the full symmetrizer of the Young operators 
for the Hecke algebra |71j . 

Applying projection operator Pa% 2 ... ap to the vectors in the tensor product Val ® (8) 

• • • ® V^ 1 , we can construct the (^ + l)-dimensional vector space Va% 2 - ae associated with the 

(2) 

spin-^/2 representation of U q (sl 2 ). For instance, we have Pa^l ^ — )a = (q/[2] ? )||2, l) a . Here 
we have introduced | H — ) a = |0) ai |l)a 2 - We denote Va^a 2 ... ae also by or Vq for short. 
Similarly, we denote Pi% 2 ... ae by Pif for short. 

(2s) (2s) 

Let us consider the tensor product Vy <g> • • • <S> , which gives the quantum space for 

(2s) 

the higher-spin transfer matrices. We construct the 6th component of the quantum space 

from the 2sth tensor product of the spin-1/2 representations: V^slb-i^+l ® " " ® ^2s(6-l)+2s' ^ or 
b = 1, 2, . . . , N s . We therefore define P^ and ?§?} L by 



7V S 7V S 

>(2a) _ TT p(2s) p(2s) _ "I T p(2s) 
12-L _ 11^2s(i-l)+l' ~~ ll- r 2s(i-l) + r 

i=l i=l 



Here we recall L = 2sN s . 
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3.6 Higher-spin monodromy matrix of type (£, (2s)® Ns ) 

Let us now introduce complete strings. For a positive integer i we call the following set of 
rapidities Xj a complete ^-string: 

A, = A - (2j - I - 1)77/2 for j = l,2,...,£. 

Here we call parameter A the string center. 

Let us now set inhomogeneous parameters Wj for j = 1,2, ...,L, as N s sets of complete 
2s-strings [TO]. We define w 2 S (l-i)+p for /3 = 1,. . . ,2s, as follows 

<4 2 s (L)+/3 = & ~ & ~ ^ for b = l,2,...,N s . (3.4) 

We now introduce the massless monodromy matrix of type (1, (2s)® s ) associated with ho- 
mogeneous grading. We define it by 

-7?(l,2s+) n ,r n ^ _ 5(2«) p(i,i+) C\ . /,„( 2s h ^p( 2s ) 

/ K 2s +)(A;{e 6 K) | {2s+) (A;{6}ivJ \ 
" ^ C^+\\;{£ b } Ns ) D(^)(A;{&U) /' 

Here, the (0,0) element is given by l( 2s +)(A; {&}atJ = P^} L A^ 1+ \X; {wf s) }l)P^). l - 

We shall now define the massless monodromy matrix of type (£, (2s)® " ) associated with 
homogeneous grading. Let us express the tensor product ® (vj 2s) ® • • • <8> V^ s) J , by the 
following symbol 

N s 



(£,(2s)® N *) = (e,2s,2s,...,2s). 

(I) (I) (I) 

Here we recall that Vq abbreviates Va 1 a a ...ai- For the auxiliary space Vq we define the massless 

monodromy matrix of type (£, (2s)® s ) by 

7?(*,2s+) _ 5W ^(l,2s+) n ^(l,2s+) / x s 

- t 0,12---AT s ~~ - r aia2---af- t ai,12---Af s V A ai^ a2j i2...AT s l A ai W " 
X -^o 2 I,12-jV s (Aai — — l)??)-Paia2-a < - 

Here we remark that it is associated with homogeneous grading. 

Let us now construct the higher-spin monodromy matrices associated with principal grading. 
From the higher-spin monodromy matrices associated with homogeneous grading we derive them 
through the inverse of the gauge transformation as follows |54| 



T {&,2sp) _ (11,2s) \ 1 T (i,2s+) (X) ( (e,2s) \ 



Here Xai—lt 12 N s denote the following: 

(1,28) _«(/) /. \*(2a)/ A N a^CA ^ 

\ „...„,. 12... v. - 9 ai ... ai (A )<P 1 (Aij ■■■<P Ns (AjvJ, 

where Ao denotes the string center, Ao = A 01 — (£ — 1)77/2. 

For an illustration, let us consider the case of £ = 1. For type (1, (2s)® Ns ) the monodromy 
matrix associated with homogeneous grading and that with principal grading are related to each 
other as follows 



th(1,2s+) , x re -, x (1,2s) rr(l,2sp) < c -, x / (1,2s) V 

r 0,12...JV a ( A ;l^}iVj = Xo,12...^ T 0,12...7V s ( A ;{^)jVj (Xo,12...ivJ 
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In terms of the operator- valued matrix elements we have 



(2a) A2sp) m f (2s) \ p -A v (2s) R (2sp) f (2s) \ 

Xl2...N s A 12-N s \ A > \X,12-N S ) e X\2-N a n 12-N a \ A l yCl2-N s ) 

X (2s) r (2sp) (x) ((2s) y 1 (2s) (2sp) m t (2s) \~ 



We shall now introduce the spin-1/2 monodromy matrices with special inhomogeneous pa- 
rameters. Let us introduce a set of 2s-strings with small deviations from the set of complete 
2s-strings 

w$$_ 1)+fl =tb-(P-l)v + er? ) for b = 1, 2, . . . , N s and = 1, 2, . . . , 2s. 

Here e is a infinitesimally small generic number and are generic parameters. We call the 
set of rapidities w^f(u-i)+p fo r @ = 1?2, ...,2s "almost complete 2s-strings". We denote by 
T'^ 1,2s+;e ' (A) the spin-1/2 monodromy matrix T^ 1 ' 1 "^ with inhomogeneous parameters wj being 
given by the set of almost complete 2s-strings: Wj = w^ s ''^ for j = 1, 2, . . . , L 

T (l,2 S +;e) m _ T (l,l+) ( X . Un {2s ' e) \ \ 

We express the elements of T^ 1,2,s+;e )(A) as follows 



ji(l,2a+;e) ^\ 



C£ £) (A) ^2 S . + £ e) (A) 



Here we recall that A\ 2 ...l (A) denotes ^4| 2 ... L (A; {u>J- ' We also remark the following: 

A?Ck (A; ) = Jim Pgi^ (A; {^-> } J pg^ . 

3.7 Series of commuting higher-spin transfer matrices 

Suppose that \£,m) for m = 0,1,. . . ,£, are the orthonormal basis vectors of and their 

dual vectors are given by {£, m\ for m = 0,1, . . . ,£. We define the trace of operator A over the 
space V® by 

i 

tr vW A= ^2{£,m\A\i,m). 

m=0 

We define the massless transfer matrix of type (£, (2s)® Ns ) by 

7{e,2s+) n ,_ (fp(t,2s+) m N 

f 12-JV s V A i - tr VW l/o,12-7V s l A J J 

= ft 1 ;^ ( A ) f i 2 1 ,itk (a -»?)•• • riX+k (a - (/ - mwtia- 



It follows from the Yang-Baxter equations that the higher-spin transfer matrices commute 

(2s) (2s) 

in the tensor product space <g> • • • <g> , which is derived by applying projection operator 
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ii2--L *° ^1 ^ ® " ■ " ® • For instance, for the massless transfer matrices, making use of ([2.8 



and (j2.9[) we show 

^IKtAl.tyM] =0 for £,m £ Z> . 

Consequently, for the massless transfer matrices, the eigenvectors of ^.jj (A) constructed by 
applying £?( 2s+ )(A) to the vacuum |0) also diagonalize the higher-spin transfer matrices, in 
particular, the spin-s massless XXZ transfer matrix, i-j^-'.AT W- Thus, we construct the ground 
state of the higher-spin XXZ Hamiltonian in terms of operators B^ 2s+ \X), which are the (0, 
l)-element of the monodromy matrix T( 1 ' 2s +) . 

3.8 Algebraic Bethe ansatz for higher-spin massless transfer matrices 

In terms of the vacuum vector |0) where all spins are up, we define functions a(A) and d(X) by 

A^(X;{ Wj } L )\Q) =a(X;{w j } L )\0), D^(X;{w,} L )\0) = d(X;{ Wj } L )\0). 
We have a(X;{wj}L) = 1 and 

L 

d(X;{wj} L ) = Y\ b(X,Wj). 

5=1 

Here b(X, /jl) = b(X — fj,). For the homogeneous grading (w = +) and the principal grading 
(w = p), it is easy to show the following relations: 

A^ 2sw \X)\0) = A^ 2sw \X)\0) = a^(X;{^ b })\0), 
D( 2sw \X)\0) = 5 (2su;) (A)|0) = d (2s \X;{£ b })\0), 

where a^^X; {&}) and d( 2s )(A;{&}) are given by 

a^(X-m) = a(X;{wf s) })=l, 



(2s)- 
,(2s)' 



d {2s) (A;U b }) = d{X;{wf s) }) = l[b 2s (X,Q. 

P =i 

Here we have defined 6j(A, fx) by b t (X, fx) = sinh(A — /i)/sinh(A — n + trj) . Here we recall b(u) = 
b\ (u) = sinh u / sinh(u + rf) . 

■ — - — (2sw) 

In the massless regime, we define the Bethe vectors |{A Q } M ) for w = + and p, and their 

— (2sw) 

dual vectors ({X a } M | for w = + and p, as follows 

In \ M (o \ ^ 
— (2sw) T ~,_ \ -- — — (2sw) r ~,„ , 

\Mm ) = U B{ 2SW) (K)\0), ({X a y M | = {0|nC (2sw) (Aa). (3.5) 

a=l ct=l 

Here we recall B^ 2s+ \X a ) = P^ s) L B( 1+ \X a , {w% ^l)?^. The Bethe vector gives an 

eigenvector of the massless transfer matrix 

for io = + and w = p with the following eigenvalue: 

A (i, 2 « 0)( j = gmh^-^ + n n smh(^ -A 3 - + ,) 

" sm^Aj-zz) y y smh^-Aj) 
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if rapidities {Xj}m satisfy the Bethe ansatz equations 



Let us denote by |{A a (e)}j^ s1i; ' e ' 1 ) the Bethe vector of M Bethe roots {Xj(e)}M for w = +,p: 

M 

7=1 

where rapidities {Xj(e)}M satisfy the Bethe ansatz equations with inhomogeneous parameters 
w^ s ' e ^ as follows 



■3 

j^ S ' e ^ as fnllnws 



q(A i (6);{ W f ;e) } L ) ^4 b(X k (e),X 3 (e)) 

It gives an eigenvector of the transfer matrix 

i^V; {^f ;e) }J = A^{»; {wf s ^} L ) + D^\^ {wf s ^} L ) 
with the following eigenvalue: 

f = i J sinh(/i- Aj(c)) 

Let us assume that in the limit of e going to 0, the set of Bethe roots {A,(e)}M ap- 
proaches {A,-}jvf- Assuming the continuity of the limiting procedure, we have 

M M 

i{A,}i; +) ) = limn (^b^kw^l) io) = pg?Lan s(3b+!e) (A,-w)io>. 

Thus, the expectation value with respect to the Bethe state of {Aj}m is given by the limit of 
that of {Xj(e)}M sending e to zero. For the B operators associated with principal grading, we 
have 

kmS* ) = an (^^HxS^j-^a^^cAy^pau^) io) 

= (xg^J^^ixi^.., x lim o nP (2sp;e) (A,( e ))|0). 

i=i 

Let us introduce symbols for the ground state of the integrable spin-s XXZ spin chain. We 
denote it by IV^ 2 ^) associated with principal grading. It is given by multiplying the projection 
operator to such a product of the spin-1/2 B operators with inhomogeneous parameters being 
given by the set of complete 2s-strings that acts on the vacuum: 

M 

\^f Sp) ) = {x^nX^IlXi^l • nP (2SKO) (A 7 )|0). 

7=1 
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We denote by \i/j^ sp '°^) the product of the spin-1/2 B operators with inhomogeneous parameters 

(2 



(2s) 

given by complete 2s-strings w - which acts on the vacuum state: 



M 

| V) (2 S p;0) ) = -Q jB (2 SP ;0) (A7) | 0) _ 
7=1 

3.9 Commutation relations with projection operators 

Let us discuss an application of the fusion construction of projection operators (|3.3|) . Hereafter 
we assume that rapidity A does not take such discrete values at which the transfer matrix 
becomes singular or non-regular, such as — V + m:\f^-Y (1 < j < L) for n £ Z [54]. Here 

we recall that w y - ' are inhomogeneous parameters forming complete 2s-strings. 

Lemma 1. Projection operators P^f and Pto?}.l commute with the matrix elements of the 
monodromy matrix ^012— { w j 2 ^}l) such as A( 2s+;0 )(A) 

pSHlI&ZI (A; {«f 'iJPgl = i-glTM.jA; {wf »},). (3.7) 
For instance we /iaue P^} l B ( - 2s +^(X)pS s } l = P§i-B (2 * +i0) (A). 

We show (|3.6p and (|3.T|) by the Yang-Baxter equation. In fact, ^^.^(A; {w^}l) com- 

(2s) 

mutes with the projection operator -P 12 ---L thanks to the fusion construction of projection ope- 
rators ([33]) [H]. We derive (J37T]) making use of OT8|) . 

3.10 Quantum inverse scattering problem (QISP) for the spin-s operators 

We can express any given spin-s local operator in terms of the spin-1/2 global operators such 
as A,B,C and D; i.e. we have the QISP formulas for the spin-s local operators [53]. For an 
illustration, we show the case of 6 = 1, i.e., we express one of the spin-s elementary matrices in 
terms of the spin-1/2 global operators. 

Lemma 2 ([53^154]). For a pair of integers i and j satisfying 1 < i, j < I, the spin-£/2 elementary 
matrix associated with principal grading is decomposed into a sum of products of the matrix 
elements of the spin-1/2 monodromy matrix as follows 

/ !\ 1/2 

{ £/3 } k=i 

Here the sum is taken over all sets of ep such that the number of integers (3 satisfying ep = 1 
and 1 < (3 < I is given by j. We take a set of e' a such that the number of integers a satisfying 
e' a = 1 and 1 < a < £ is given by i. We have expressed the element of (a,/3) in the monodromy 
matrix T^ e )(X) by T^J p;e \x) for a, (3 = 0, 1. 

For an illustration, let us consider the spin-1/2 formula [H [11] (see also |10t I73|): 

n— 1 n 

= II • II (3-8) 

i=i i=i 
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Here we recall that the spin-1/2 transfer matrix t^^ mL (X) is given by the trace of the monodromy 

matrix of type (1, 1® L ): t$}. L (X) = A^ lp \\) + D^(X). We remark that the expression M 
holds if inhomogeneous parameters Wj (1 < j < L) take generic values. Multiplying the expres- 
sions of formula (|3.8I) for n = 1, 2, . . . , m, we have 



• • • = tr (ej L (u;i))tro( e S ,J ' a 4^... L (^)) • • • 

m 

x tr (e^<i.LKn)) II ( 3 - 9 ) 

3=1 

Here we note that we have i?Q^f (0) = ^-o,n from the normalization condition of the R- matrices, 
where Ho, n denotes the permutation operator acting on the 0th and nth sites (see also Sec- 
tion [33]). Thus, we have 

3=1 



We note that the QISP formulas (|3.8[) hold if the inhomogeneous parameters are generic. If 

(2s) 

we send them to a set of complete 2s-strings such as w - , then the transfer matrix becomes 
non-regular or singular, and relations such as ()3.10j) do not hold. Instead of complete 2s-strings, 
we therefore put "almost complete 2s-strings", w[ 2s ' e \ into inhomogeneous parameters Wj. Here 

3 J 

parameters uf? s ' € ' are generic, and hence the QISP formula (I3.8D holds. 

3.11 Expectation value of a local operator through the limit: e — > 

In the massless regime, we define the expectation value of product of operators E^'^ 2 

k=l 

-(2sp) 

with respect to an eigenstate |{A a } M ) by 

ft ^-«|{X>r> 

(II Et h(2 '») ({A.(g"") = Jl% lp) _ (2ap , ■ (3.11) 

k=1 ({^a} M \{^a} M ) 

In order to evaluate (|3.1ip we make use of the following formulas. 

Proposition 2 ([53\ I54j). Let us take a pair of integers i\ and ji satisfying 1 < < £. For 
arbitrary parameters {^o}n and {\/3}m with i\ — j\ = N — M we have 

N M 

(oi n ■ e\ imip) ■ n ^ p) (A^io) 

a=l 0=1 



\ 



ji 



-IN M 

X>| n C^°\» a ) ■ ef' E1 ■ ■ ■ ef et ■ [J B^(\p)\0). (3.12) 

1 {ep} a=l 0=1 

Here we take the sum over all sets of £p such that the number of integers (3 with er = 1 for 
1 < f3 < £ is given by ji . We take a set of e' a such that the number of integers a satisfying 
e' a = 1 for 1 < a < I is given by i\. Each summand is symmetric with respect to exchange of e' a ; 
i.e., the following expression is independent of any permutation tt £ S^: 

N M 

<0| n C^ ;0) (/ia) • ep 1 • • • ep< • J] B^°\\p)\0). (3.13) 

a=l 0=1 



One-Point Functions of the Integrable Spin-1 XXZ Chain 



19 



Here we remark that S n denotes the symmetric group of n elements. 

We evaluate the expectation value of a given spin-s local operator for a Bethe-ansatz eigen- 

-~-(2s) 

state ({Aq,}^ ), as follows. We first express the spin-s local operators in terms of the spin-1/2 
local operators via formula (|2.10p . Through Proposition [2] the expectation value of the spin-s 
local operators is reduced into those of the spin-1/2 local operators. We now assume that the 
Bethe roots {X a ( e )}M are continuous with respect to small parameter e. It follows from the 
assumption that each entry of the Bethe eigenstate |{Afc(e)}j^ s,<E ' ) ) is continuous with respect 
to e. Then, we apply the spin-1/2 QISP formula with generic inhomogeneous parameters w^ s '^ 
such as formula (|3.9p . We next calculate the scalar product for the Bethe state |{Afc(e)}j5 S,< ^)- It 
has the same inhomogeneous parameters w^ s,e ^ as the global operators appearing in the QISP 
formula, so that we can make use of Slavnov's formula of scalar products for the spin-1/2 case. 
Calculating explicitly the determinant of the scalar product with Slavnov's formula, we can 
show that the expression of the scalar product is continuous with respect to e at e = 0. Thus, 
sending e to 0, we obtain the expectation value of the spin-s local operator (|3.11j) . 

Corollary 1. Suppose that i\ and j\ are integers satisfying 1 < < £, and {fJ.k}N a- re 

arbitrary parameters. Let us assume that Bethe roots {A 7 (e)}jv/ are continuous at e = and 
lim A 7 (e) = A 7 for 7 = 1,2,..., M. We have the following: 



N M 
k=l 7=1 











9 





-1 
g 



x 



N M 

J2\™(0\UC^X,, k ).T^^ ■ n^(A 7 ( £ ))|0). (3.14) 

{e } € k=l 7=1 



Here we take the sum over all sets of e^s such that the number of integers (3 satisfying eg = 1 
for 1 < f3 < t is given by j\ . We take a set of e' a such that the number of integers a satisfying 
e' a = 1 for 1 < a < t is given by i\. We have defined 4> m ({X 7 }) by <^ m ({A 7 }A/; {vjj}i) = 

m M 

J b(X a — Wj) with b(u) = sinh(u)/ sinh(u + n). 

j=l a=l 

We can evaluate the form factors and the expectation values of a spin-£/2 operator through 
Corollary Q] [53] . The corrections of the form factors given in the paper |19] are listed in referen- 
ce [20] of the paper [54] (see also [53]). Corrections for the paper [20] are listed in reference [21] 
of the paper [54] . 

For an illustration, let us consider the spin-1 case. We calculate the one-point function 
(El' 1 ^}. Here we have i± = j\ = 1. Setting s[ = and e' 2 = 1, we have 

(4 2P) I^' 1(2P V? P) ) = (^ 2p;0) |ere 2 'Vf ;0) ) + (^le^e^V^ 05 )- 

Here we have taken the sum over sets {ep} such as {ei = 0,^2 = 1} and {e\ = 1,£2 = 0}. 
Making use of the spin-1/2 QISP formula we have 



pp;e) , (2;e) 
'12-L\ W j 
3=1 

e o,i e i,o = c^iwf^B^iw^) n (t?£i(wf'' e) 
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Therefore we have 



(E 1 ) - 2 ({A 7 }, {«,, } J hm 

(4 2p;e) | C(^) (wf e) )B(^ (wf e) )\ 4 2p;e) ) \ 



+ lim 



~° (4 2p;e) i4 2p;e) ) / 

4 Quantum group symmetry relations in the spin 1 case 

We show some important topics. We derive symmetry relations among the expectation values 
of products of the spin-1/2 operators from the spin inversion symmetry. In particular, we show 
how to transform the basis vectors constructed in the 2sth tensor product space of the spin-1/2 
representations to the 2s + 1-dimensional vectors in C 2s+1 . 

4.1 Rotation symmetry of the XXX spin chain and irreducible components 
of operators 

Let us consider the XXX case where the SU(2) symmetry holds for the total spin operators. 
The tensor product of two spin-1/2 representations of si (2) decomposes into the direct sum 
of spin-1 and spin-0 representations; i.e., V(l) (g> V(l) = V(2) V(0). Here we recall that 
V(l) <8> V(l) is four-dimensional, and the components V{2) and V(0) are three-dimensional and 
one-dimensional, respectively. In the spin-1 representation V(2) we have the basis vectors and 
basis covectors as follows 

112,0) = || + +), (2, 0|| = {+ + ||, 

|| 2 ,i> = || + -> + !!-+>, (2,i|| = i((+-|| + (- + ||), 

112,2) = || --), (2,2|| = (--||. 

Here we recall that || — h) denotes |l)i ® 1 0) 2 - 

In the spin-0 representation V(0) we have the basis vectors and basis covectors as follows 

||o,o) = 1 + ->-|-+}, (0, o|| = !((+- || -(- + 11). 

In terms of the basis of the spin-1 irreducible representation we express the symmetric projection 
operator as follows 

P {2) = ||2,0)(2,0|| +||2,1)(2,1|| + ||2,2)(2,2||. 

In the spin-s XXX case we define elementary matrices by 

E m,n(2s) = ||2 a , m )(2a,n||. 

In the tensor product V(l) £3 V(l) there are sixteen elementary matrices e^ ,kl e^ ,k2 for 
jl,j2,ki,k2 = +• For an illustration we express the operator e\ + in terms of the basis 
vectors and their covectors as follows 

(2,l|| e +-e 2 - + ||2,l) = i (0,0|| e +^e2 + ||2,l) = i, 

(2,l|| e +-e2 + ||0,0) =-i (0,0||e+-e2 + ||0,0) = ~. 
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In terms of the bases of vectors and covectors, we have 

ef-e^ + = i(||2,l)(2,l|| + ||0, 0){2, 1|| - ||2,1)(0,0|| - ||0,0)<0,0||). (4.1) 
Applying the projection operators to the right-hand-side of (|4. 1 1) . we have 

P( 2 )e+-e 2 -+p( 2 ) = i||2,l)(2,l|| = l -E^ 2 \ 
Similarly, we have 

e++e 2 — = i(||2,l)(2,l|| + ||0,0>(2,l|| + ||2,l)(0,0|| + ||0, 0)(0, 0||) . 
We thus have 

P^ef+e^-p^ = i||2, 1){2, 1|| = -E l ' 1{ ~ 2 +\ 
In terms of irreducible components, we have 

P( 2 )e+-e 2 - + p( 2 ) = pWeret + P (2) - 
We thus have 

pWe^ + e+-P {2) = P^ef-e^ + P^ = P^ef+e^P^ = P^eret + P {2) - 

We shall evaluate the expectation values of spin-s local operators by reducing them into those 
of the spin-1/2 local operators. Applying formula (|2.10p to the case of e' x = and e' 2 = 1, which 
correspond to + and — , respectively, we have 



<^ 2) I^ 1(2 W> = <^ 2;0) l^ 2 - + l^ 0) > + <4 2 ;°>| e + + e2 -|#°)>. 

Here we remark that the vector \ijjg ) is given 

by|4 2S;0) >= I! 5 (2s;0) (A 7 )|0), while the vector 

7=1 

\ipg 2s ^} is given by multiplying the projection operator: = P\2 l n B^ 2s '°\\y)\0) . 

7=1 

4.2 Spin inversion symmetry 

For even L we may assume the spin inversion symmetry: U\ip g 2sp ' ^} = ±|t/^ 2sp ' ^) for U = f\ a j- 

3=1 

Here we recall that associated with the ground state of the integrable spin-s XXZ spin chain the 
vector |4 2SP;0) ) is given by |4 2SP;0) > = U £ {2sp;0) (A 7 )|0>. 

7=1 

We derive symmetry relations as follows |53[ [54] 

{4 2sm \4' £l ■ ■ ■ 4r e2s \4 2sP '' 0) ) = (4 2sp;0) |e;^ ,1_£l • • • e 1 -^' 1 -^^ 2 ^). (4.2) 

Applying the spin-inversion symmetry (j4.2j) we derive symmetry relations among the expec- 
tation values of local or global operators [53l [54] ■ 

For an illustration, let us evaluate the one-point function in the spin-1 case with i\ = j\ = 1, 
^i,i(2p)^ Setting e\ = and e 2 = 1 we decompose the spin-1 elementary matrix in terms of 
a sum of products of the spin-1/2 ones 

(4^\eI^\4^) = (4^%r4M 2p '' 0) ) + rt 2p;0) l^ 1 4'V? p;0) )- 
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Through the symmetry relations f)3. 13|) with respect to e' a we have the following equalities: 

<4 2p;0 ^°4'Vf ;0) > = (^ 2p;0) |e^ e°'V? ;0) >, 
{^%Yef\^) = (^ ;0) |er4'°l4 2p;0) )- 
From spin-inversion symmetry (|4.2p we have 



<^ p; V4'Vf ;0)) > = (Vf^l^'Vf ;0) >, 
(^ »|e? 1 4' |^°)) = <V^ 0) l4'°4'Vf ;0) > 
and hence we have the equalities of the four terms. We therefore obtain the following: 

{^f)\E\^f)) = 2<^|er4'Vf ;0) >- 

We thus derive the double-integral representation of the one-point function (El' 1 ^ 2 ^) of [20], as 
we shall show in Section [6l 

4.3 Transformation from V^ 2 ^ to the (2s + l)-dimensional vector space C 2s+1 

We shall express the spin-s massless XXZ transfer matrix (2s + 1)^ x (2s + 1) N ° matrix 
acting on the tensor product of the (2s + l)-dimensional vector spaces C 2s+1 ; i.e., acting on 

In Section 13.61 we have defined the spin-s XXZ transfer matrix through the fusion method. 
It is expressed in terms of operators defined on the Lth tensor product space of the spin-1/2 
representations, (y^')® L , and given by a 2 L x 2 L matrix. We have constructed them by applying 
the projection operators to the spin-1/2 XXZ transfer matrix with inhomogeneous parameters 
given by complete strings Wj . 

We now formulate the spin-s XXZ transfer matrix in terms of the basis of the (2s + 1)- 
dimensional vector space C 2s+1 such as |2s, to)) for m = 0, 1, . . . , 2s. As the basis vectors of the 
(2s + l)-dimensional representation of U q (sl2) we introduce vectors \2s,m) with the following 
normalization: 




2s, m) = ||2s,m)/W I I for m = 0, 1, . . . , 2s. 



We denote by (2s, m\ the transposition of |2s,m) 

(2s, m\ = (|2s, m)) 1 for m = 0, 1, . . . , 2s. 

Let us denote the complex conjugate of a complex number z by z. We express the Hermitian 
conjugate of a vector |2s,m) by 

(2s, m\ = (|2s, m)y for m = 0, 1, . . . , 2s. 

Let us introduce the transformation 5: y( 2s ) — y C 2s+1 . We define it by 

2s 

S = |2s,m))(2sT 

m=0 

It maps the basis of the spin-s representation V^ 2 ^ constructed in the tensor product space 
® ■ ■ ■ (g) V^'i i-e., 1 2s, m) for to = 0, 1, ... , 2s, to that of the (2s + l)-dimensional represen- 
tation C 2s+1 ; i.e., |2s, to)) for m = 0, 1, . . . , 2s. We can show the following relations: 

ggi,i(2sp)s^ = for *,j =0,1,..., 2s. (4.3) 
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Here we recall that E 1 ' 3 denote the (2s + l)-by-(2s + l) unit matrices which have only one nonzero 
element 1 at the entry of (i, j) for i,j = 0, 1, . . . , 2s. 

For an illustration, let us consider the spin-1 case. The basis vectors of C 3 are given by 

|2,0)) = (1,0,0)*, 12,1)) = (0,1, 0)*, |2,2)) = (0,0,1)*. 

In the spin-1 case the transformation S: — > C 3 is given by 



S=£|2,m))<2, 



■/?? 



m=0 



In the massless regime where q is complex with \q\ = 1, explicitly we have 
10 

T2 * 

1 
Taking the Hermitian conjugate of S we have 



S 



St = £|2,ro>((2, 



/ 1 \ 

o -4 o 



m 



m=0 











V2q 

V o o i / 



It is straightforward to show the following: 



1 
SS f = | 1 o I , 
1 

In terms of the bras and kets we have 

2 2 



f 1 








\ 





1 

2 


1 

2 








1 

2<2 


1 

2 





^° 








1 ) 



SS*= ^|2,m))(2,m|^|2,n)((2,n| 

m=0 n=0 
2 2 2 

= ^2S(m,n)\2,m))((2,n\ = ^ |2,m))((2,m|. 



m=0 n=0 

Similarly, we have 

2 



m=0 



5 f 5= ^|2,m)(2, 



771=0 



In order to transform the conjugate vectors ||2,m) it is also useful to introduce the complex 
conjugates of transformations S and S^: 



S=X>,m))(2, 



771=0 



10 
1 
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5 f 



|2,m)((2,m| 

m=0 



/ 1 





\ 

75 ° 



V o o i / 

They are related to the projection operator P^ 2 \ We have 

/ 1 \ 

I JL n 2 

^ %m)(2,m\ 



| 

§ 
\ 



j_ 

2</ 

1 

2 





p(2)_ 



m=0 



The spin-1 elementary matrices E l '^ 2p ^ are transformed into the 3x3 unit matrices E % ^ as 
-ggi,i{2p)-g\ = E i,j for ^ j = o, 1, 2. 

For instance we have 







£gl,l(2p)3t 



S 1 = £ 



i.i 




We have thus confirmed relations (14.31). 



Let us introduce the transformation which maps the tensor product of the spin-s represen- 



tations: V, 



(2s) 



® • • • ® vj^ s ^ to the tensor product of the (2s + l)-dimensional representations: 
Hi w 

Si ® ■ ■ ■ ® : V^' 
We also define its complex conjugate 



q2s+1 



C A f + . We define it by the tensor product of transformation S as follows 



r (2s) 



q2s+1 



Si ® ■ ■ ■ ® 5tv 3 : ^ (2s) ® ■ ■ ■ ® -> C? 



2s+l 



<-i2s+l 



Let us consider the spin-s ground state with (2s + l)-dimensional entries, |\I/g S ^). For the spin-1 
case, it gives the ground state of the spin-1 XXZ Hamiltonian (jl.l|) . In terms of the ground 



state constructed by the fusion method, \ip. 



[2sp\ 



it is given by 



\* { G S) )=Sl® 



A2sp) s 



Here we recall that \ipg ) denotes the ground state of the integrable spin-s XXZ spin chain 
constructed through the fusion method, where the evaluation representations are associated with 
principal grading. In terms of the eigenvector with (2s + l)-dimensional entries, the expectation 
value of a given local operator E with (2s + l)-dimensional entries is given by 

<*g s) |£|*g s) ) = (Vf P) \S\ ® • • • ® S^ESi ® • • • ® S Na \4 2s ri). 

Therefore, the operator E corresponds to the operator E^ 2sp ^ in the fusion construction as follows 

E^ = s\ ® • • • ® S^ESi ® ■ • • ® 5jv s - 
For instance, from (|4.3p we have the following: 



jgi,i{2sp) _ g^^Jg 



for i,j = 0, 1, 



2s. 



Similarly, we have the following relations for the spin-s XXZ transfer matrices defined as (2s + 



1)^ x (2s + 1)^ matrices t 



12-JV.' 



to those of the fusion construction: 



7(£,2sp) 
L 12-N a 



(A) 



31 



sk*?2*l(A)sl® 



® SV 



for 



1,2, 



2s. 
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5 Multiple-integral representations for spin-s case 

We introduce some useful symbols for expressing the correlation functions of the integrable 
spin-s XXZ spin chain. We derive the multiple-integral representation of the spin-s correlation 
functions by mainly following the procedures of [20] except for the formula of reducing the 
higher-spin form factors into the spin- 1/2 scalar products such as in Corollary [TJ 

Let us sketch the main procedures for deriving the multiple-integral representation of the spin- 
s XXZ correlation functions. First, we introduce the spin-s elementary operators as the basic 
blocks for constructing the local operators of the integrable spin-s XXZ spin chain. Secondly, 
we reduce them into a sum of products of the spin-1/2 elementary operators, which we express 
through the spin-1/2 QISP formula in terms of the matrix elements of the spin-1/2 monodromy 
matrix, and evaluate their scalar products through Slavnov's formula of the Bethe-ansatz scalar 
products. Here, the expectation value of a physical quantity is expressed as a sum of the ratios 
of the Bethe-ansatz scalar products to the norm of the Bethe-ansatz eigenvector. Furthermore, 
the ratios are expressed in terms of the determinants of some matrices. Thirdly, solving the 
integral equations for the matrices in the thermodynamic limit, we derive the multiple-integral 
representation of the correlation functions. 

Let us summarize the multiple-integral representations of correlation functions for the inte- 
grable spin-s XXZ spin chain in a region of the massless regime with < £ < 7r/2s |20j. We 
show the revised expression [53\ 154] . Here we recall that in the massless regime we set rj = 
with < C < 7T. 

We express any correlation function of the integrable spin-s XXZ chain in terms of the linear 
combination of products of the spin-s elementary matrices. They are defined by 

m 

Fi 2sp H^,j k }) = (4 2s *\ n^ Jfc(2sp Vf p) )/(4 2sp) l4 2sp) )- 

k=l 

where E % ^^ k<y2sp ^ denotes the (2s + 1) x (2s + 1) elementary matrix whose entries are all zero 

except for the (ik,jk) element which is given by 1, for each k with 1 < k < m. Here integers ik 

and jk satisfy 1 < ik,jk < 2s. We recall that IV^ 2 ^) denotes the spin-s ground state. 

Let us consider a product of the spin-£/2 elementary matrices associated with principal 

grading, ^ ■ ■ ■ £^ I '"' m ^ p \ for which we shall evaluate the zero-temperature spin-s XXZ 

\kV \k\ 

correlation functions. We introduce variables e a and which take only two values or 1 for 
A; = 1,2, ... ,m and a, f3 = 0, 1, . . . , I. For the mth product of elementary matrices, we introduce 
sets of e^a s and e^'s (1 < k < m) such that the number of as satisfying = 1 and 1 < a < 2s 

\k] 

is given by ik and the number of /3s satisfying = 1 and 1 < /3 < 2s by jk, respectively. We 
then express them by integers e'jS and EjS for j = 1, 2, . . . , 2sm as follows: 

= £ a Y for a = 1,2,... ,2s; fc = l,2,...,m, 

\k] 

e 2s (k-i)+p = £p for /3 = 1,2,... ,2s; k = 1,2,... ,m. 

We express the mth product of (2s + 1) x (2s + 1) elementary matrices in terms of a sum of 
2smth products of the 2x2 elementary matrices with entries {ej,e^}; i.e., e^ 1 ' 61 ■ ■ ■ e £ 2 2 s s ™' £2sni 
[201 [54]. 

For given sets of Ej and e'j for j = 1,2,..., 2sm we define a~ by the set of integers j satisfying 
e'j = 1 and a + by the set of integers j satisfying Ej = 0: 
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We denote by a_ and a+ the number of elements of the set a and a + , respectively. Due to 
the "charge conservation", we have 

a_ + a + = 2sm. (5-1) 

m m 

Precisely, we have a_ = ^ i\. and a+ = 2sm — ^2 jk- Here we recall that for the .R-matrix of 

fc=i fc=i 
the XXZ spin chain matrix elements R(u)^ vanish if a + b ^ c + d, which we call the charge con- 
servation. It follows from the charge conservation that the correlation function ({sj, e' }) 

m m 

vanishes unless the two sums are equal: ^2 ik = 3k- We therefore obtain relation (|5.ip . 

k=i k=i 

We remark that the charge conservation of the i?-matrix corresponds to the "ice rule" of the 
six-vertex model, which is defined as a two-dimensional ferro-electric lattice model. 

For sets a~ and a + we define A-,- for j £ a~ and X'j, for f £ a + , by the following sequence: 

( A k ax ' • • • ' A Ln' Ajmin ' ' ' ' ' Xj ™^ = ( A i' • • • ' A 2sm)- 

Let us recall the assumption that in the region < £ < tt/2s the spin-s ground state IV^ 2 ^) 
is given by N s /2 sets of the 2s-strings: 

X[ a 1 = fia -( a -l/2)r ] + 5i a) for a = 1, 2, . . . , N s /2 and a = 1, 2, . . . , 2s. 

Here we also assume that string deviations 5^ are very small for large N s . In terms of rapidities 
forming strings, A« , the spin-s ground state associated with the principal grading is given by 

N 3 /2 2s 

iv? sp) > = II II B(2sp) ( A i Q) 'ia}7vj|o). 

a=l a=l 

Here we have M Bethe roots with M = 2sN s /2 = sN s . The density of string centers, p(X), is 
given by 

P(A) 



2Ccosh(7rA/C) 



which has simple poles at A = i£(n + 1/2), for n € Z with the residues (— l) n /(27ri). 
We define the (j,k) element of a matrix S = S[(Xj)2 S m'i (^j^^sm) by 

Sj,k = p{Xj-wf s) +r l /2)5(a{X j ),P(k)) for j,k = 1,2, ... , 2sm. 

Here 5(a, f3) denotes the Kronecker delta. We define (3(j) by 

m=3-28[[(j-l)/2s]], 

where the Gauss symbol [[x]] is defined by the greatest integer less than or equal to a real 
number x. We define ct(Xj) by a(Xj) = 7 (1 < 7 < 2s) if Xj is related to the integral variable 
fij by Xj = [ij — (7 — 1/2)7/, or Xj « where j3{k) = 7 (1 < 7 < 2s) [2D]. We remark that /x,- 
correspond to the centers of complete 2s-strings Xj. When we evaluate a(Xj), we assume that 
the integral paths of J^_^_i^± ie are replaced by those of J^'^ly^ for 7 = 1, 2, ... , 2s, 
respectively. Here we remark that when we deform the integral path of /f^^J ^sjjf i6 to that of 

X^oo-i(7-i/2X ^ or ^ = ^' ^' ■ ■ ■ ' we ma y bave the contribution of a simple pole at A = w^ 2 ^ 
with integer k satisfying f3(k) = 7. 
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With the above notations, we express correlation functions for the massless spin-s XXZ chain 
in the form of multiple integrals as follows 

(poo+ie poo — if+ie poo — i(2s— l)£+ie \ 

/ +/ +•••+/ dAi-. 

J— oo+ie J— oo — if+ie J — oo — i(2s— l)f+ie / 

oo+ie poo— i(2s— l)f +ie \ 

+ •••+/ dA Q+ 

-oo+ie J — oo—i{2s—l)C+ie I 

oo— ie p'OO — if — ie poo— i(2s— l)f— ie \ 

+ / +•••+/ dA Q++ i • • • 

-oo— ie V— oo— if— ie J — oo— i(2s— l)f — ie / 

oo— ie poo—i(2s—l)^—ie\ 

H h / dA 2sm 

-oo— ie J — oo— i(2s— l)f — ie / 

x ^2 Q({ e 3' £ 'jh^T ■ ■ ,A 2sm )det5(Ai,. . . ,A 2sm ). (5.2) 

Here we have defined Q({ej, e'j}; Ai, . . . , A 2sm ) by 
Q({ £ ji £ 'j}'i -Vi, - ■ ■ , A 2sm )) 

n nsinh(A i -4 2s) +r ? ) n sinh^- 



(_!)<*+ j£Q- \fc=l fc=J+l 



n sinh(A^ — A fc + 77 + e^ fc ) 

l<fc<-£<2sm 

n sinh(A;,-4 2s) -r/) n sinh(A;,-4 2s) ; 
fc=i fc=j'+i 

sinh — 

Kk<K2sm 



(5.3) 



Here we have set e ki e = ie for Im(Afc — A^) > and e^e = — ie for Im(A/ c — A^) < 0, where e is an 
infinitesimally small positive real number: 0<e<l. The normalization factor C is given by 



c({i k ,j k }) = n {^F(i,i k )/F{i,j k ) q ^-^i 2 ^^i 2 ) = n 

6=1 \ 



k=l 



where q = e 71 = e 1 ^. 

Here we should remark that in (|5.2p the sum of a + ({Ej}) is taken over all sets {ej} corre- 
sponding to {ejg'} (1 < k < m) such that the number of integers f3 satisfying ej^ = 1 with 
1 < /3 < 2s is given by j^, for each satisfying 1 < A; < m. In (|5.3|) we take a set a~({e'j}) 

corresponding to for fe = 1, 2, . . . , m, where the number of integers a satisfying = 1 and 
1 < a < 2s is given by i k for each k (1 < k < m). 



6 Multiple integrals of the spin-1 one-point functions (s = 1) 

We calculate analytically the integrals for the spin-1 one-point functions. Considering the 
residues which are derived when we shift the integral paths, we explicitly evaluate the dou- 
ble integrals expressing the spin-1 one-point functions. Hereafter, we shall often denote the 
spin-s elementary matrices E l '^ 2sp ^ by E 1 ^ for simplicity. 
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6.1 (E 22 ): The emptiness formation probability 

Let us evaluate the emptiness formation probability (EFP) (E 2 < 2 ( 2 p)}. In this case we have 

k=h = 2; ( £l ,e 2 ) = (l,l), (4,4) = (1,1); C = l; 
a + = 0, a" = {1,2}; (Ai, A 2 ) = (Ai, A 2 ). 

Here the symbol denotes the empty set. We evaluate EFP as follows 



(^ 2(2p) ) = 2 ({A 7 };{^}J(lim 



(2)" 



(4 2p;e) iD( 2 p^)( w f ;£) )£>( 2 ^ £ )(4 2;£) )i4 2p;£) ) 
<4 2p;e) i4 2p;e) > 



Let us denote the integral path J^^^ ia by f c . The multiple-integral formula reads 



22\ 



+ 

C_ ie VC. 



dA 2 Q(Ai,A 2 )det5'(Ai,A 2 ), 



where Q(Ai, A 2 ) and 5(Ai, A 2 ) are expressed in terms of (p(x) = sinh(x) as 
9?(Ai - ^ 2) )(^(A 2 - ^i 2) + 7]) 



,( 2 V 



<5(Ai, A 2 ) — 

ip(X 2 - Ai + 7/ + e 2 i)(/2(u4 - u> 2 ) 

5(Ai A 2 ) = | p ^ Al ~ w * + n ^ 'Wi)- 1 p ( Xl ~ ™ 2 + V/ 2 ) ^(Ai),2 I 
\p{\ 2 - wf> + rj/2) 5 a (A 2 ),i p(A 2 - ?4 + i]/2) (5 a ( A2 ) i2 / 

We now shift the integral paths C-\ t and C- v -i e into C_ v / 2 and C_ 3j? / 2 , respectively. During 
the contour deformation each of the integral paths does not cross any pole of the integrand, and 
hence we have 

(E 22 ) = I [ + f ) dAi ( / + / ] dA 2 Q(A 1 ,A 2 )det5(Ai,A 2 ). 

We now denote C_ v / 2 and C_ 3r) / 2 by C± and C 2 , respectively. After expanding the above 
expression with respect to the types of integral paths, we have four terms. However, only two 
of them survive due to the Kronecker deltas in the matrix S 



(E 



22 \ 



p(X 1 - wf ] + 7?/2) 
p(A 2 - w {2) + T//2) 

p{X 1 - w {2) + 77/2) 

p(A 2 - 4 2) + ri/2) 

p(Ai - w {2) + 77/2) 

p(A 2 - w {2) + 77/2) 

p(Ai - 4 2) + ??/ 2 ) 



dAx / dA 2 Q(A!,A 2 ) 



+ / dA x / dA 2 Q(Ai,A 2 ) 
JCi Jc 2 

+ f dXi f dA 2 Q(Ai,A 2 ) 

JC 2 JCi 








+ / dAi / dA 2 Q(Ai,A 2 
Jc 2 JCi 



p(A 2 - w {2) + 77/2) 

dAi / dA 2 Q(Ai, A 2 )p(Ai - w {2) + v/2)p(X 2 - w [2) + t?/2) 



Ci JC 2 

^ dAi I dA 2 Q(Ai, A 2 )p(Ai - 4 2) + W2)p(A 2 - w? 5 + 77/2), 
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(21 (2) 

Substituting w\ = £i, w 2 = £i — r), we have 



(E 22 ) 



d//i / d/i 2 (Qi2 - <22i)pOi - ii)p{^2 - ii) 

J — oo 
/>oo 

dxi / dx 2 (Qi2 - <52i)p(a;i)p(a;2) 



where x\ = p\ — £i, x 2 = /^2 — £i 5 and Q12 and Q21 are given by 

1 ^(/xi - £1 + 77/2)^2 - 6 - W2) 



Ql2 = Q(w-??/2,/X 2 - 377/2) 

1 <^(xi + r]/2)p(x 2 - ry/2) 



^(/x 2 - Mi - ie) 



¥>(??) ^(^2 - si - ie) 



(6.1) 



Q21 = <9(mi - 377/2, /j 2 - 77/2) 



1 ^ - ei - 77/2)^(7x2 - a + W2) 



(^(77) </?(/j 2 - pi + 2r/ + ie) 

1 (f(xi - T]/2)if(x2 + 7//2) 



<p(v) ^2 ~ xi + 2t/ + ie) 
Thus, we have 

(E 22 ) = J12 — /21 
where Ji 2 and I21 are given by 



I12 
hi 



dxi 



-OO 
OO 



dxi 



dX2Ql2(xi,X2)p(x 1 )p(x 2 ), 

dx 2 Q 2 i(xi, x 2 )p(x 1 )p(x 2 ). 



The integrand Q12 is transformed into Q21 when we shift the integral path as x\ — > x\ — r] 
and x 2 — > X2 + 77. First we shift the integral path in /12 as x\ — > x\ — 77. 



Xx 2 -ie 

-n/2 

fj VWVWWW WWWWWWWWWWWV\A>- 

Here we note that due to the sign in front of e in the denominator of (|6.ip . the integrand Qi 2 
has a pole at x\ = x 2 — ie as a function of x\. Here we recall that e is an infinitesimally small 
positive real number. We therefore express the integral I12 in terms of a sum of two integrals, 
J 1 + J2, as follows 



'12 



dx 2 p(x 2 ) I dx 1 Qi 2 (x 1 ,x 2 )p(x 1 ) 

D 

OO — TJ 
OO — r; 

dx 2 p(x 2 ) I dxxQuixi - rj,X2)(-l)p(xi) 



dx 2 p(x 2 ) 



dx 1 Q 12 (x 1 ,x 2 )p(xi) - 27riRes [Qi 2 (xi, x 2 )p(xi)\ xi=X2 _ ie ] 



2tt\ I dx2p(x 2 )Res[Q 1 2(xi,x 2 )p(x 1 )\ Xi __ 



Ji + J 2 - 
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Here we have made use of the anti-periodicity: p(x + nrj) = (— l) n p(x). We also remark that 
the simple pole at x\ = —rj/2 due to p{x\) is canceled by the factor y?(xi + r]/2) in Qi 2 (xi,X2). 

Let us first consider the single integral J 2 derived from the pole at x\ = x 2 — e. Explicitly 
evaluating the integral J2 we have 



/oo 
dx 2 p(x 2 
-00 



2 



1 ip(x 2 + r)/2)ip(x 2 - rj/2) 

00 <p{v) -1 

7r f°° , cosh 2x — cos C C — C cos C 
dx- 



4C 2 sin C Loo cosh 2 (7rx/C) 2( sin 2 C 

Here we have made use of formula (|B.2p . 

\X2 

X>7/2 



Let us next consider the double integral J\. We shift the integral path in J\ as x 2 — > x 2 + rj. 
We derive the wanted integral l 2 \ as follows 



00 i-oo 



Ji = / dxi / dx 2 Qi 2 (xi - rj,x 2 )(-l)p(xi)p(x 2 ) 

-oo J — OO 

00 roo+r) 

dxi / dx 2 <5i2(xi - r/,x 2 )(-l)p(xi)p(x 2 ) 

-OO J —OO+T] 

/OO /'OO 
dxi / dx 2 Qi2(xi - 7?,x 2 + r/)(-l) 2 / o(xi) y o(x 2 ) 
-00 J — 00 

/oo /*oo 
dxi / dx 2 Q 2 i(xi,x 2 )p(xi)p(x 2 ) = I 2 i. 
-00 J —00 

Here we note that the simple pole at X2 = 77/2 due to p(x 2 ) is canceled by the factor (p(x 2 — rj/2) 
in Qi 2 . 

Finally we have the analytical expression of the one-point function (E 22 ) as follows 

/P 22v_, T _ (T , n r _ 7 _ C-srnCcosC 
; — J12 — l 2 i — [Ji + J 2 ) — l 2 i — J 2 — — _ . 9 — • 

2£ sin C 

6.2 = 2( e i' 1 e°'°) 

Let us calculate a spin-1 one-point function, (E 11 ). Setting = 1 and e 2 = in formula ()2.10|) 
we have 

n=Ji = l; (ei,e 2 ) = (0,1), (1,0); C = 1; 

a+ = {l}, {2}; a"={l}; (Ai, Ai) = (A l5 A 2 ), (A 2 , Ai) = (Ai, A 2 ). 

Here (ei,e 2 ) = (1,0) corresponds to a + = {2}, and hence we have (A' 2 , Ai) = (Ai, A 2 ). 
Applying formula (|3.14p we express as follows 

(^ 1,1(2p) ) = (^i^^^Vrv^i^) 

^ 2 ( { A 7 },K }^hm iiUe^Ue)}^) 



One-Point Functions of the Integrable Spin-1 XXZ Chain 



31 



({A a ( e )}^ e) pfe^)(4 2;e ))Afe^)(4 2;e ))|{A a ( e )}^ e V 



lim " ^ /Jj " 1 ^ / IL " v ' JJM ' . (6.2) 



Considering the spin inversion symmetry and the quantum group invariance we evaluate (E 11 ) 
by 

(^ 1 ' 1(2P) ) = 2^( { A 7 };{^} L ) 



({A a ( e )}if e) |Z?^-)(4 2 - ) )A(^.)(4 2 -))|{A Q ( e )}^ ;e) ) 



x lim \i-'u\-jim v~i >_ \-2 n^±X21M ' (6 3) 

<{Ac(e)}S P;e) |{Ac(e)}if e) > 



Let us briefly review how we reduce (j6.2|) to (|6.3p . It follows from formula ()3. 12[) that we 
have 

(4 2p) i^' 1(2p) i4 2p) > = (^iwi^ ) + (^ki^ri^ ). 

Due to the spin inversion symmetry (|4.2[) we have 

We have the following symmetry relation (|3.13|) due to the quantum group invariance £/g(sZ 2 ): 

{^W** ;0) ) = {^le^e^Vf ;0) }- 
Therefore, we have 

(^IWI^ ) = (^ 2p;0) l e i' le 2' l4 2p;0 >- 

Here we also recall that \^g P '^) = |{Aa(e)}^ 2p '^)- 

Let us consider the case of (si,s' 2 ) = (1,0) and (ei,e 2 ) = (1,0), which corresponds to a + = 
{2} and ex = {1}, and hence we have (A 2 , Ai) = (Ai,A2). The multiple-integral formula reads 

(E U ) = 2U +1 )dAi(/ +f ]dA 2 Q(A 1 ,A 2 )det5(A 1 ,A 2 ), 

where Q(Xi, A 2 ) and det 5(Ai, A 2 ) are given by 

n(x , v y(A 2 - ^ 2) )y(Ai - wf ] - rj) - - r/)<p(A 2 - ^ + 77) 

V?(A 2 - Ai + 7] + e 2 i)pK 2) - u4 ') <P{M - A 2 - 77 + £12)^(77) 

5(Ai A 2 ) = I P ^ Al ~ ^ + ^ <5 "( Al )' 1 P( - Al ~ ^ + v ^ "Wo.* ) 

\p{\ 2 ~ wf ] + 77/2) <5 a (A 2 ),l P(A 2 - ?4 + 77/2) <5 a (A 2 ),2 / 

Here we recall that Jq denotes the integral path f^^"* a and also that ip(x) = sinh(x). 

ia (2) 
Let Tj denote a small contour rotating counterclockwise around A = w - for each j. We 

shift the integral paths C-\ e — > C%, C- v -\ e — > C 2 and C + \ t — > C\ — F\, C- v+ i t — > C 2 — T 2 , where 

C\ = C_ ri / 2 and C 2 = C_ 3r; / 2 . For instance, we have 



L 



/ dAi. 



Expanding the determinant of matrix S, we thus obtain 



(E n )/2 =( [ - [)d\ 1 [ dA 2 Q(Ai, A 2 )p(Ai - 4 2) + ^/2)p(A 2 - 4 2) + W2) 
KJd Jr 1 J Jc 2 
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C*2 JT% 



dAi / dA 2 Q(Ai,A 2 )p(Ai -w { 2 2) +r]/2)p{\ 2 - wf ] +77/2). 



The one-point function (E n ) is now expressed in terms of J±, J 2 , i£i and K 2 , as follows. 

(^ 11 ) = 2(-^ 1 + K 2 + J 1 - J 2 ). 

Here we shall give definitions of integrals J±, J 2 , K\ and K 2 and calculate them shortly in the 
following. For K\ and K 2 , making use of the formula: 27riRes [p(X — w + rj/2)\ \ =w ] = 1, we 
have 

Kt = f dAi f dA 2 Q(Ai, A 2 )p(Ai - + 7?/2)p(A 2 - w { 2 2) + v /2) 

JT-l JCi 

/oo 
dp 2 Q(^,p 2 -3rj/2)p(p 2 -Ci) 
-00 

dxp(x), 



C 2 

00 



and 



K 2 



d\ l / dA 2 Q(Ai, A 2 )p(Ai - w { 2 2) + r?/2)p(A 2 - w { ? + rj/2) 



Ci 



/ dA 2 Q(4 2) , A 2 )p(A 2 - tt^ + r?/2) = / dp 2 Q(Ci - r,, p 2 - r,/2)p(p 2 - 

J Ci J~oo 



2 cosh 77 / dxp(x) 



<p(x + 7]/2) 



2 cosh 77 / dxp(x 



ipjx - Tj/2) 
ip(x + 77/2) 



<p(x + 3r?/2) 
We have defined the integrals J\ and J 2 by 

Jx = I dAi / dA 2 Q(A 1> A 2 )p(Ai - w {2) + r,/2)p{\ 2 - + rj/2) 

JCi JCo. 

1 \ tp(x 1 - 3rj/2)ip(x 2 - rj/2) 



Ci 
00 



Jo 



c 2 



dxi / dx 2 p(xi)p(x 2 ) 

^(t?) / ip{x 2 -x 1 -\e) 
dA a / dA 2 Q(Ai, A 2 )p(Ai - + ??/2)p(A2 - wf ] + r?/2) 

1 \ <p( Xl - 5r]/2)<p(x 2 + r?/2) 



Ci 



dxi / dx 2/ o(xi)/9(x 2 ) 



vC 7 ?) / ( / J ( a; 2 - xi + 2r/ + ie) 



As in the case of (E 22 ), we transform the integral J\ into J 2 by shifting the integral path as 
x\ — > X\ — rj and x 2 — > x 2 + rj. First we shift the integral path in J\ as X\ — > x\ — rj. There 



are two simple poles at x\ = x 2 — ie and x\ = —rj/2. Using 27riRes p(x)\ x= _ r) , 2 
calculate the residues as 



T, we can 



27riRes 



p(xi)p(x 2 
-2(cosh r])p(x 2 ) 



1 \ ip(xi - 3rj/2)<p(x 2 - v/2) 



and 



2-7riRes 



p(xi)p(x 2 ) 



<p(rj) J <f(x 2 -x\- ie) 

<f(x 2 - rj/2) 
(p(x 2 + 77/2) ' 

1 \ tp(xi - 3ri/2)tp(x 2 - rj/2) 



Xl=-Tj/2 



<p(v) 



ip(x 2 — x\ — ie) 



xi=X2— ie 
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2vri 

<p{v) 



p(x 2 ) 2 ip(x 2 - 3r]/2)ip(x 2 - rj/2). 



Thus we have 

Ji = ~h ~h + 



OO — 7? 



dxi 



h + 



-oo—r] 

oo 

dxi 



where 



h 
h 



dx 



-oo 



-2(coshr])p(x 



dx 2 p(xi)p(x 2 ) 

oo 

dx 2 {-l)p{x 1 )p{x 2 ) 



ip( X - 7]/2) 



1 \ (p(xx - 3rj/2)(p(x 2 - 77/2) 



<p(rj)J tp(x 2 -x 1 -ie) 

1 \ ip(xi - 5r]/2)(p(x 2 - t?/2) 



<p(rj) J <p(x2 - xx + rj - ie) 



(p(x + rj/2) 



dx 



2vri 

Mv) 



p(x) 2 f(x — 37]/2)tf(x — rj/2) 



Next we shift the integral path as x 2 — > x 2 + rj. Here we remark that the simple pole at x 2 = rj/2 
of p(x 2 ) has zero residue due to the factor (p(x 2 — rj/2) of the integrand. Thus we have 



-h -I 2 + 
h -I2 + J2 



dx^ 



dx 2 (-l) 2 p(xi)p{x 2 ) 



1 \ ip(xi - 5ri/2)ip(x 2 + rj/2) 



<p(rj)J ip(x 2 - x x + 2rj - ie) 



where we have omitted the infinitesimal e since we can shift the integral path without crossing 
the poles. Thus, we have 

(E u ) = 2{-Kx +K 2 + Jx- J 2 ) = 2{-Kx + K 2 - h - I 2 ) = 2{-Kx - h), 

where we have used the fact that K 2 = I\. Using the formula (jB.ip . we have K\ = 1/2. Next 
we consider the integral I 2 . Shifting the integral path of x as x — > x + m, we have 



oo+ii" 



tp(x — rj/2)tp{x — 3^/2) 



-00 + 17T 

+ 27riRes 



p(x) 2 dx 



<p(x - y/2)ip(x - 377/2) 2 

v(v) P[x) 



x=m/2 



ip(x + rj/2Mx-rj/2) 2 2 ip(-rj) 
<p( V ) ( l) P{X)dX + ^\(rj)(2^ 

sinh(x + rj/2) sinh(x — rj/2) 



4CV(?7) 



COsh 2 (7Tx/C) 



-dx 



2vri 



Making use of the formula: sinh(x + rj/2) sinh(x — rj/2) = (cosh2x — cosh7y)/2 we have 



2vri 



cosh 2x 
cosh 2 (7rx/C) 



dx — cosh 77 



1 



cosh (irx/C,) 



-dx 



C — sin C cos ( 
~ 2Csin 2 C ' 



where we have used the formula (|B.2|) . Finally, we obtain 

cos C(sin C — C cos C) 



(E u ) = 2(-Kx-I 2 



Csin £ 
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6.3 (E 00 ) 

In this case we have 

h=3i = 0; (ei,e 2 ) = (0,0), e' 2 ) = (0, 0); C=l; 

a + = {l,2}; a~ = 0; (A' 2 , Ai) = (Ai, A 2 ). 

The multiple-integral formula reads 



(E 



00 \ 



/ + 



dAi | / +/ dA 2 g(A 1 ,A 2 )det5(A 1 ,A 2 ), 



where Q(Ai,A 2 ) and S"(Ai,A 2 ) are given by 

V?(A 2 - u4 2) )¥>(Ai - iuf - ??) 



Q(Ai,A 2 ) 



S(X 



(2) (2) ' 

ip(\ 2 - Xi + rj + €2i)<p(w\ - w 2 ) 

7J 

(2) (21 

^?(A 2 - + r//2) 5 a (A 2 ),i P(A 2 - w\ } + rj/2) 5 Q ,(A 2 ),2 / 



, /p(Ai -w™ +r//2)<5 a(Al)i i p(Ai - u»J^ + T//2) 6 a{Xl) 2 

1, A 2 J — / ^ (2 ) i ,„ ' ' ("21 . 



(2) 



Here we recall that J* c . denotes the integral path X^ !'j Q and V 9 ( x ) = sinh(x). 

We now shift the integral paths C + \ t — > C\ — Y\, C- v+ \ e — > C 2 — T 2 , where C\ = C_ ri / 2 , 

(2) 

C 2 = C_3 v /2 an d Tj is a small contour rotating counterclockwise around A = Wj . Expanding 
the determinant of matrix S, we obtain 



^ \lci iri) 



I ~ I ) dA 2 Q(Ai, X 2 )p(X 1 - w {2) + r]/2)p(X 2 - w { 2 2) + V /2) 

C 2 JTu 



dAi I / - / 1 dA 2 Q(Ai, A 2 )p(Ai - + ?7/2)p(A 2 - + rj/2) 
'c 2 \JCi ./rv 

Jj- J 2 _/g_ J 4 + / 5+ / 6 



where 



/!=/" dAi / dA 2 Q(Ai,A 2 )pip|, / 2 = / dAx/ dA 2 Q(A 1 ,A 2 )^p? ) 

I 3 = Q{w { 2 2 \w {2) ), h= f dX 2 Q{w {2 \x 2 )p 2 2 , h= [ dX 2 Q{w { 2 2) ,X 2 )pl 

JC 2 JCi 

h = I dA 1 Q(A 1 ,^ 2) ) /0 2' 

and pi = p(Xj — w k + t?/2). Shifting the integral path as for the former cases, we have 
1 1 — I 2 = —K\ — K 2 where 



A", 



dx 2 27riRes 



1 tp(x 1 -3r}/2)ip(x 2 -ri/2) 
cp(rf) y?(x 2 -xi- ie) 



p{xi)p{x 2 ) 



X1=X2—1£ 



sin £ cos C — C 
2C sin 2 C 



+ 1, 



K 2 



dx 2 27riRes 



1 (p(xi - 3r]/2)ip(x 2 - 77/2) 
99(77) if(x 2 - xi- ie) 

ip(x - 7]/2) 



p{x{)p{x 2 ) 



xi=-'n/2 



2 cosh 77 / dx^ — _^.=ip( z ). 



One-Point Functions of the Integrable Spin-1 XXZ Chain 



35 



The other terms are calculated as 

/oo roo 
p(x)dx, I 5 = K 2 , h = ~ J p(x)dx. 
-oo J — oo 

Summing up all the contributions, we have 

<o = c ;; inC 2 7 c . (6.5) 

2Q sin Q 

Here we can confirm that the relation (E 22 ) = (E 00 ) by directly evaluating the integral. 
6.4 (E 11 ) through 2(e° 1 ' 1 e^°) 

We evaluate (E 11 ) by calculating the multiple integral representing (e^e^' ). Here we recall 
that due to the spin inversion symmetry we have (E 11 ) = 2{e^ 1 e^)- In this case we have 

n=ji = l; (ei,e 2 ) = (0,1), (e^) = (1,0); C = 1; 
a+ = {2}; a ={2}; (A' 2 , A 2 ) = (Ai, A 2 ). 

The multiple-integral formula reads 

(E U ) = 2U +1 ]dAi(/ + / ]dA 2 Q(A 1 ,A 2 )det5(Ai,A 2 ), 

where Q(Xi, A 2 ) and det 5(Ai, A 2 ) are given by 

q(\ 1 a 2 ) = ^( A 2 - w i 2) + ^M^i - wT ~ v) _ <f(*i ~ 6 - ??)y(A 2 - 6 + ??) 

(p(X 2 - Ai + 7] + e 2 i)^(u;f ) - w { ?) <^ Al ~ A2 ~ 11 + ^M??) 

5(A X A 2 ) = | P ^ Al ~ ^ + ^'WlP ~ *4 + W 2 ) *a(Ai),2 ] 

\p(A 2 -wf ] +v/2)8 a (\ 3 ),i P( A 2-^ 2 2) +W 2 )^(A 2 ),2/ 

Here we remark that we have the same Q(Ai,A 2 ) as in equation (I6.4p for the case of (e^e^' )- 
We therefore obtain 



11 cos£(sin£ — £cosC) 



C sin 2 £ 



We have thus confirmed the quantum group invariance (e^e?,' ) = (e}' 1 ^' ) through the 
multiple- integral representation. 

Finally in Section [6] we give an important remark: through an explicit evaluation of the 
multiple integrals of (E 1 ' 1 ^ 2 ^) we have shown the following relations: 

/ 0,0 1,1a / 1,1 0,0\ / 1,0 0,1\ / 0,0 1,1\ 

(e{ e 2 ' ) = (e{ e 2 ' ), (e{ e 2 ' } = (e{ e 2 ' ). 

It follows that in the spin-1 case, every one-point function is expressed in terms of a single 
multiple integral, which corresponds to the expectation value of a single product of the local 
spin-1/2 operators. In general, however, the spin-s correlation function of an arbitrary entry is 
expressed in terms of the expectation values of a sum of products of the local spin-1/2 operators 
such as shown in (|5.2p . Here we recall that the sum over sets a + ({es}) in (|5.2|) corresponds to 
the sum over sequences {eg} in the reduction formula of Corollary [TJ 
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<E 22 > = <E 00 > a 




/ 

<E U > 

i 

Figure 1. One-point functions of spin-1 XXZ chain in the massless regime obtained by the explicit 
evaluations of multiple integrals. The red and blue lines represent those for (E 22 ) = (E 00 ) and (E 11 ), 
respectively. 



7 Consistency with numerical estimates 
of the spin-1 one-point functions 

We now show that the analytical expressions of the spin-1 one-point functions are consistent with 
their numerical estimates, which are obtained by the method of numerical exact diagonalization 
of the integrable spin-1 XXZ Hamiltonian. 

Let us fist summarize the analytical results derived in Section [6l Evaluating the multiple 
integrals explicitly, we have obtained all the one-point function for the integrable spin-1 XXZ 
chain as 

/g2,2(2p)\ = /go,o(2 P h = C ~ sinCcosC /_gi,l(2p)\ = cosC(sinC- CcosQ 

2£ sin 2 ( ' £ sin 2 £ 

which are shown in Fig. [TJ In particular, via evaluation of the multiple integrals, we have 
confirmed the uniaxial symmetry relation: 

(E 22 ) = (E 00 ). (7.1) 

Through the direct evaluation of the multiple integrals we confirm the identity: (E 22 ) + + 
(E 00 ) = 1. Here we recall that assuming the uniaxial symmetry (|7.ip the analytical expression 
of (E 00 ) has been given in |20| . 

Furthermore, we have confirmed the relations among the correlation functions due to the 
quantum group U q (sl2) symmetry and the spin inversion symmetry as follows 
= 2^4,1) = 2(e id e o >0) = 2(e o,i e i, 0) = ^1,0^ 

In the XXX limit A -> 1 we have (E 22 ) = {E n ) = (E 00 ) = 1/3, which has been shown by 
Kitanine in the XXX case [IT]. In the free Fermion limit A ->■ we have (E 22 ) = (E 00 ) = 1/2, 
and (E 11 ) = 0. Here we should remark that we consider the region < £ < 7r/(2s) with s = 1, 
namely, < A < 1. 

Finally, we confirm the analytical results by comparing them with the numerical results of 
exact diagonalization, which are shown in Fig. [2] In Fig. [2] the red and blue lines represent 
the analytical results obtained by evaluating the multiple integrals of the one-point functions, 
(E 22 ) = (E 00 ) and (E n ), respectively. The black dotted lines represent the numerical estimates 
of the one-point functions which are obtained by the method of exact diagonalization of the 
integrable spin-1 XXZ Hamiltonian with the system size of N s = 8. We numerically obtain the 
ground-state eigenvector of the integrable spin-1 XXZ Hamiltonian, and calculate the numerical 
estimates of the one-point functions. 

We have found that the numerical and analytical results of the spin-1 one-point functions 
agree quite well in the region < A < 1, as shown in Fig. [2j We thus conclude that the 
numerical results should support the validity of the multiple-integral representations for the 
spin-1 one-point functions. 
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Figure 2. Comparison with the exact numerical diagonalization. The red and blue lines represent 
analytical results obtained by the multiple integrals for (E 22 ) = (E 00 ) and (E 11 ), respectively. The black 
dotted lines represent those obtained by exact diagonalization with the system size N s = 8. 



A Derivation of reduction formula (12.101) 



For the spin-£/2 Hermitian elementary matrices associated with homogeneous grading, E i, ^^ ,+ \ 
we introduce coefficients gij by 



Then, we have 



'£ 




'£ 




i 


q 


j_ 





J(i-l)/2-j(j-l)/2 n -(a(l)+---+a(i)-i)+(b(l)+...+b(J)-j) 



We derive the reduction formula for the Hermitian elementary operators £ li >-W+) as follows 



'£ 




'£ 


."'( 


i 


q 


J_ 




a a(i) 


\£,0)q- 



'£ 



J(i-l)/2-i0-l)/2 



IM 



(A.l) 



{zp}e {e' a }e 

Here {£p}i is given by a sequence of or 1 such that the number of integers j3 for 1 < f3 < £ 
satisfying eg = 1 is given by j, and {e' a }e is given by a sequence of or 1 such that the number 
of integers a satisfying e' a = 1 is given by i. 
We can show the following property: 

Lemma 3 ([M])- Let ot be a set of distinct integers {a(l) , . . . , a(i)} satisfying 1 < a(l) < 
• • • < a(i) < £, we have the following: 



i\\o~„ 



-(a(l)H ha(i))+i 



n -1 



i(i-\)/2 



'a(l)--- (T a(^ ) ( l 

which is independent of the set a~ = {a(l),a(2), . . . ,a(i)}. 

Applying Lemma [3] we show that the inside of the parentheses (or the round brackets) of 
equation (jA.ip is independent of a(k)s. Making use of the following: 



£ 1 
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where e' a are such a sequence of or 1 that the number of e' a = 1 is given by i. We thus have 

-l 

£*d'(4+) 







'£ 




i 




J_ 






q 







^-l)/2-,0--l)/2|| M(M | 



_t(i-l)/2-i(j'-l)/2| 



Here we recall that {£/3}i is a sequence such that the number of integers j3 of 1 < j3 < £ 
satisfying ep = 1 is given by j. The integers a(k) (1 < A; < i) and (1 < k < j) satisfying 
1 < a(l) < • • • < a(i) < I and 1 < 6(1) < • • • < b(j) < £, respectively, are related to the 
sequences {e' a }e and {e^i by the following relation 



£,,£! 10 



1,0 0,0 0,0 0,1 



„o,i 



e a(l) "' e a{i) e l "' e l e b(l) "'%■)• 



B Useful integral formulas 

r0 ° dx 



oo coshx 
00 cosh2ax 



7T, 



dx 



cosh 2 x 
Acknowledgments 



2-7ra 



sm 7ra 



for lal < 1. 
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